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ABSTRACT 

An analysis of the measurement accuracy requirements of a high- 
resolution meteor radar for observing short -period atmospheric waves 
is presented, and a system which satisfies the requirements is described. 
A medium-scale, real-time computer is programmed to perform all echo- 
recognition and coordinate-measurement functions. The measurement 
algorithms are exercised on noisy data generated by a program which 
simulates the hardware system, in order to find the effects of noise 
on the measurement accuracies. Under typical conditions, the height 
accuracy is ±1 km, and the velocity accuracy ±5 msec these accuracies 
are sufficient for the study of the structure of atmospheric winds in 
the 80-120 km height region. 
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1. INTRODUCTION 

1.1 Winds in the Upper Atmosphere 

The neutral winds in the atmosphere can be decomposed into three 

components: prevailing wind, tidal wind, and irregular wind [Beer, 1972]. 

The prevailing wind is predominantly east-west. This east-west flow varies 

through the course of the year in a characteristic pattern in both hemispheres 

according to season. Diurnal periodic motions due to atmospheric tides are 

superimposed on these prevailing winds. Tides are most visible in wind 

studies at 85-105 km [ Greenhou and Neufetdj 1956]. Tidal theory predicts that 

- 1/2 

the amplitude of these waves should increase upwards as p , where the 

neutral gas density p decreases almost exponentially in the atmosphere. 

This increase in amplitude is a direct consequence of conservation of energy. 

If A is the amplitude of an atmopsheric tidal wave, then the kinetic energy 
2 

is given by p A / 2 , which must remain constant in the absence of dissipative 
forces. 

Since the atmosphere is stratified by a gravitational, field, there 
exists a class of waves in which the air particle motions are transverse 
relative to the direction of propagation of the wave. When the effects 
of gravity and density stratification become dominant, atmospheric gravity 
waves result. These atmospheric gravity waves are so named because their 
behavior is influenced by gravitational forces. Such waves have interested 
meteorologists for quite a long time [Hauruitz, 1951] . Observation of winds 
reveals a scattering of velocities about the regular oscillations of the 
tidal wind. This irregular wind makes the largest contribution to the 
total wind vector, and increases in amplitude with height, just as tidal 
waves do. Hines [1960] showed that these apparently random variations could 


be explained in terms of internal gravity waves. Greenhcw and Neufeld 
[1959] found that the irregular winds in their 1956 data had a time scale 
of about 100 minutes and vertical scale of about 6 km. Photographs of long- 
enduring meteor trails provide a more direct measurement of the vertical 
scale of the irregular winds. These trails are rapidly distorted from 
their initial straight-line shape. In recent years, most of the information 
about vertical structure of the irregular winds has come from chemical 
trails released by rocket experiments. Observation of these trails indicate 
unambiguously that the irregular winds are predominantly horizontal, with 
amplitudes increasing with height. 

Much of the irregular motion in the ionosphere was at one time thought 
to be turbulence. Turbulence power estimates (i.e., the rate of dissipation 
of turbulence energy) were as high as 25 watts per kilogram [ Booker , 1956] ; 
such high estimates were based on motions now attributed to gravity waves. 
Turbulence is much weaker than was previously believed, but is nonetheless 
present at heights up to 100 km [Hines , 1963] . 

Turbulence is obvious from the diffusive growth of meteor trails. 

The trail cross section increases first because of molecular diffusion, 
but in a few seconds eddy diffusion gains importance and ultimately 

dominates. Greenhao [1959] has deduced from this transitional growth 

-3 

that a turbulent power of 7 x 10 watts per kilogram is present at 
meteor heights . 

With this value for turbulent power, Hines [1963] attempted to get 
some idea of the scale of the turbulent spectrum [Booker , 1956] . Values 
for the length scale L^, time scale T and velocity V ^ of the small scale 
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end of the spectrum can be calculated once the kinematic viscosity factor 

2 -1 

is known. This factor is approximately 10 m sec at 95 km [A iinzner 3 et al. 3 
1959]. The values were calculated, yielding 

L 2 = 20 id 
T n - 40 sec 
V 2 = .5 msec" 1 

This scale of L n is comparable to the scale observed by Blamont and Jager 
[1961], and is also compatible with the scale of ionization irregularities 
that are observed by the scattering of radio waves \Bowtes 3 1959] . It is 
important to note that, in fact, the scale of turbulence is not a constant. 
The corresponding values for the large-scale end of the spectrum (X^, 

7 1 ) are not uniquely specified. They are, however, limited by the processes 
which generate the turbulence. These processes are generally believed 
to be associated with the irregular wind. Hines [1963] chose < 1000 sec, 
and then calculated V ^ < 3msec 1 and L ^ < 3000 m. 

These various estimates are based on the application of standard tur- 
bulence formulas. Since these formulas are derived for incompressible 
fluids, and ignore anisotropies present in a stratified fluid such as the 
ionosphere, they can be taken as conditionally acceptable relationships which 
may or may not be verified through further study. 

1 . 2 Measurement Techniques 

Wind measurements above 80 km fall into two general categories; those 
in which the experimenter places a probe or tracer in the atmosphere; and 
those in which the experimenter relies on some natural phenomenon to act 



4 


as a tracer. In the first category, several techniques are available, among 
which are grenades, Pitot tubes, falling spheres, and chemical releases. 

Each of these techniques uses a fairly large rocket to carry the experiment 
to the required height, and each returns usable data for a rather short 
period, ranging from a few minutes to an hour. All of these rocket 
techniques are expensive, and limited by the short measurement interval. 

Several natural phenomena are useful for atmospheric wind measurement. 
Observation of noctilucent clouds is one of the oldest wind measurement 
techniques [Witt, 1962], Noctilucent cloud occurrence is unpredictable, 
however, and limited to twilight hours. For these reasons, the method 
is not useful for systematic wind studies. Mitra [1949] is responsible for 
a technique which tracks the movement of ionospheric irregularities by 
ground-based radio measurement. This method is, however, subject to 
difficulty in interpretation of fading records. 

Observation of meteor trails is another wind determination method. 
Manning et. at. [1950] showed that meteor trails move with the surrounding 
atmosphere, and are therefore, useful in wind studies. Visible meteor 
trails occur rather infrequently, but ground-based radar systems can 
detect the invisible but much more frequent trails left by meteors smaller 
than a grain of sand. Since meteors enter the atmosphere throughout the 
day, a nearly continuous sampling of atmospheric wind is available from the 
method . 

Ionospheric motions have also been observed using partial reflection 
systems with spaced receivers. The velocity of a diffraction pattern formed 
on the ground by radio waves partially reflected from the 60-100 km region 



is measured by comparing signals received at separated antennas. Meteor 
radar, however, is a more common measurement technique, and the technique 
with the most potential for continuous wind study. 

1.2.1 Meteor-radar technique. A meteor entering the atmosphere is 
heated by frictional forces. The high temperature of the meteor causes 
ionization of the atmosphere it is passing through, and also causes self- 
vaporization. The ionized molecules from the atmosphere and the meter 
itself form a trail in the wake of the meteor. This ionized trail is a 
good specular reflector for VHF radio waves; that is, VHF waves are 
reflected from the trails, with the angle of incidence equalling the 
angle of reflection. A radio wave is reflected back toward its source if 
the meteor trail is oriented perpendicular to the direction of propagation 
of the wave. 

Meteor trails move with the surrounding atmosphere, causing a Doppler 
frequency shift in the reflected wave. For a pulse transmitted at a 
frequency reflected from a target moving with velocity v relative to the 

Is 

radio source, the received signal / is given by the expression 




1 + v/q 


1 - v/c J t 


( 1 . 1 ) 


where c is the velocity of propagation of the radio wave. For the condition 
V much less than o , the received signal may be approximated 
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Wind velocities in the atmosphere are certainly much less than the speed of 
light, so that the approximation is valid. The Doppler frequency shift 
of the returned signal is 

A C 1 - 3 ) 

Once formed, meteor trails grow in size through diffusion. As the 
trails grow, reflected signal initially increases in strength and then 
decreases to zero as the trail is dissipated. 

The majority of meteor- radar stations measure nothing more than the 
Doppler shifts of meteor trail echoes. Since meteor trails occur in the 
height region from 80 - 120 km, such systems measure an average wind. This 
average is weighted, however, because most returns come from about 95 km. 

More elaborate systems measure the angles of arrival and range of the re- 
flection point in addition to the Doppler shift. With this information, the 
height of the reflection can be calculated, and a wind profile of the 
80 - 120 km region can be developed. 

Other schemes for measuring height have been used. One technique is 
to observe the rate of decay of the received signal amplitude. Since the 
atmospheric density decreases with height, meteor trails at higher heights 
diffuse faster than those at lower heights. This method can yield reflection 
heights to within ±2 km at best [ Lee and Geller 1973] . 

1.2„2 Limitations of previous systems. Existing meteor -radar 
experiments have a number of limitations. Because of uncertainties in 
height measurement, wind shears and sharp wind velocity gradients are 
smoothed by spatial averaging, and cannot be accurately measured. Also, 
long delays occur between the collection of raw data and the development 
of wind profiles. These delays are due to the fact that collection occurs 
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at one place, and processing at another; in addition, processing has 
traditionally involved many manual operations. For these reasons, experi- 
mental programs which use the results of one experiment to plan another 
are severely limited. 

The rate of meteor trail echoes has also been a limiting factor. 

The observation of gravity waves requires that hundreds of meteor trails 
per hour be recorded, while current systems can provide only tens of trails 
per hour [ Barnes , 1968]. Another requirement for gravity wave observation 
is height accuracy greater than 1 km. Such accuracy has not yet been 
achieved. 

1,2.3 Advantages of the -proposed system. The Illinois meteor-radar 
system has a number of advantages over other systems. Height uncertainty 
is reduced to ±1 km by improvements in elevation angle and range measurement 
schemes. A medium-scale computer is included as an integral part of the 
system so that the velocities and polar coordinates of specular reflection 
points can be determined in real time. Much of the measurement and control 
hardware has been replaced by programmed algorithms, bringing improvements 
in overall system reliability and flexibility. Lastly, the system 
transmitter has a power output of 5 Mw, which is about two orders of 
magnitude more power than other meteor -radar systems used. This high 
power will yield about one order of magnitude more meteor trail echoes 
per unit time than other systems. 

1.3 Statement of the Problem 

The research reported on in this thesis covers four areas: 

1) Specification of free parameters such as transmitter 

pulse repetition frequency (PRF) , pulse width and shape, and 


receiver bandwidth. 
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2) Design of a software system for the collection and processing 
of data, and allocation of the computer attributes to the tasks 
involved. 

3) Design of the required measurement algorithms. 

4) Demonstration of the data reduction routines through computer 
simulation of the hardware system and collection software. 

Design decisions made in any one of the first three areas has a direct 
impact on the other two. The role of the fourth area of effort is to 
analyze the total effect of any decision relating to the first three areas. 
The over-all effort of this research, therefore, is to recognize the 
tradeoffs involved in each design decision, and maximize the usable data 
output of the system, subject to the accuracy and hardware constraints. 

1.4 Outline of Thesis 

The remainder of this thesis is organized as follows. Section 2 
describes the atmospheric phenomena to be observed, and gives the system 
constraints required to measure the phenomena, such as measurement 
accuracy and required echo rate. Also, the constraints set by the system 
which were pre-defined. These system features are the given, around which 
other aspects must be planned. Section 3.2 is the beginning of the 
description of the system design effort. Sections 3.3 and 3.4 detail 
the software data organization, and outline the levels of processing 
involved. The logical flow of software control is also given. Section 4 
describes each of the required routines. In addition, the effect of 
quantization and sky noise on algorithm accuracy, and the impact of 
hardware parameter selection (e.g. PRF, pulse width) on data reduction, 
is examined. Section 5 presents the system capabilities and signal-to-noise 
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ratio requirements. Section 6 gives a summary of all work done, and presents 
some ideas for future efforts. In particular, a re-evaluation of the 
software-hardware tradeoffs is made, and a proposal for the addition 
of some external control hardware to the system is given. 



2. SCIENTIFIC CRITERIA AND SPECIFICATIONS 


2 . 1 Scientific Objectives 

The objective of the Illinois meteor-radar program is to build a 
system that is compatible with nationally and internationally coordinated 
meteor wind programs, and also with specific research efforts where a single 
station is of value. Some of the observational programs requiring a 
coordinated network are investigation of planetary wave propagation into 
the lower thermosphere, investigation of winter changes in mesospheric 
and lower thermospheric circulation, observational studies of the latitudinal 
structure of tides, and studies of the geographic variation in atmospheric 
turbulence. Observational programs of scientific interest for a single 
station include studies of diurnal and seasonal variation of turbulence, 
comparison of thermospheric circulation and solar activity, and observation 
of the small-scale wind structure. 

Figures 2.1 and 2.2 [Rosenberg , 1968] show the north-south and east-west 
components of a collection of 60 wind profiles measured by observation of 
chemical releases from rockets. Rosenberg resolved the data shown in 
these two figures into total wind speed shown in Figure 2.3, and wind shear, 
shown in Figure 2.4. The most that could possibly be expected of a meteor 
radar system is that it produce measured wind profiles comparable to those 
from rocket releases. 

Rosenberg and Zimmerman [1972] studied the shear-related properties of 
winds in the meteor region. Table 2.1 contains some of their results. From 
these data, measurement criteria for the meteor-radar system were developed. 
It is desirable to build a system that can: 
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Figure 2.1 Figure 2.2 


Wind component velocity versus altitude from [Rosenberg, 1968]. 
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Table 2.1 

Shear -related properties of winds in the 90-120 km region 


1. Altitude, km 


90-97 97-103 103-110 


2. Total samples 


101 


534 


764 


3. Mean velocity, ms 


-1 


48 


58 


60 


4. Wave amplitude, ms 1 35 


45 


48 


5. Shear,. , ms * km 1 
lim 


75 


75 


75 


6. Shear , ms * km 1 
avg' 


21 


21 


21 


7. Shear , ms" 1 km 1 
rms 


27 


27 


27 


8. Wavelength , km 14 

6 mean 


16 


18 


9. Wavelengthj^ m> km 


3.5 


4.0 


4.5 


110-120 


482 


64 


52 


50 


14 


18 


24 


6.0 



IS 


1) Measure winds in excess of 100 msec \ with accuracy a small fraction, 
of typical atmospheric wave amplitudes. 

2) Resolve shears on the order of 50 msec * km \ 

3) Measure height, with accuracy a small fraction of the vertical 
dimension of the atmospheric waves of interest. 

There are other requirements for the longer period dynamics. These 
include measurement of two horizontal components, and the capability for 
extended operation. 

2.2. Number of Meteors Per Hour 

The shortest period oscillation that can be observed by a meteor 
radar is determined by the meteor echo rate. The Brunt-Vaisala frequency, 
u> , is the highest frequency wave that can exist in an isothermal atmosphere. 
This frequency is given by the relation 

■»> = (y _ i) 1/2 g/C (2.1) 

y 

c 

Y = = 1-4 

y 

-2 

g = 9.8 msec 
C = speed of sound 

Figure 2.5 adapted from Georges [1967] shows that the period (2tt/oj ) 

Q 

of the Brunt-Vaisala frequency in the 80-120 km region is about 5 min. 
Observation of this frequency requires at least one return per min from 
any given height. If meteor echoes were uniformly distributed over the 
80-120 km region, observation of this wave structure with one km height 
accuracy would require 40 echoes per min, or 2400 per hr. Meteor echoes 
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Figure 2.5. Period of the Brundt-Vaisala as a function of 

height, calculated using the 1962 U, S. Standard 
Atmosphere (adapted from Georges [1967]), 
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are not uniformly distributed, however, as is shown in Figure 2.6 
[Barnes, 1972], A meteor echo is most likely to occur at 98 km; 60% 
of all echoes occur at 95-99 km, and 90% at 90-100 km. Although these 
exact figures vary for different meteor-radar stations, the general 
shape of the distribution is maintained. Based on this distribution, 
the following conclusion is made: About 1000 echoes per hr are needed 

in order to resolve short-period waves near 95 km. Meteor-radar study 
of such short -period waves over the entire 80-120 km region requires 
a much higher echo return rate, on the order of 10,000 per hr. Lower 
frequency oscillations could be discriminated over a large height range 
without this high return rate, however. 

2.3. Measurement Accuracy 

The study of short-period gravity waves requires that the height, 
azimuth, and velocity of specular reflection points on meteor-trails be 
measured within certain tolerances, as discussed in Section 2.1. Since 
meteor-radar stations measure polar coordinates, the tolerance requirements 
for height and horizontal velocity measurements must be restated in terms 
of range, elevation angle, azimuth angle and radial velocity accuracies. 

In this section, each of the required system accuracies is considered in 
turn, and requirements for the polar coordinate measurement accuracies are 
developed. 

2.3.1 Height accuracy requirements. Neglecting earth curvature, the 
height of a specular reflection point on a meteor trail is given by the 
relation 

h = R sins (2.2) 

where R is the radial range, and s the elevation angle (Figure 2.7). From 
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Figure 2.6. Meteor count versus altitude at Stanford [Barmes , 1972] 
where U is for east-west and V is for north-south winds. 



TRANSMITTER 


Figure 2.7. Illustration of the measurement coordinates. 
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Section 2.2, the most likely height of a meteor-trail reflection is about 
98 km. The antennas for the Illinois meteor radar are aimed with their 
main lobes elevated at 45°. The representative radial range for meteor- trail 
echoes is therefore 98/sin 45°, or about 140 km. 

The principal contribution to the total energy reflected from a 
straight-line meteor trail comes from the first Fresnel zone. This zone 
is the interval on the trail between the point of specular reflection, and 
the points where the range has increased by X/4 over that to the specular 
point, where A is the radar wavelength. The normal range differential of 
\/2 for the end points of a Fresnel zone is halved in this case because 
of the round-trip path of the radar pulse. The width of this zone is 
v0j?/2 measured from the specular point along the trail. Because of 
symmetry, the zone from which the radar beam is reflected can be taken 
as a segment on the meteor trail of length centered on the specular 

point. For the Illinois meteor-wind system, the radar wavelength is 
7.33 m. Using 140 km as the range, a representative dimension for the 
first Fresnel zone is calculated to be about 1.4 km; the vertical pro- 
jection of this zone is about 1 km. The Doppler shift of a meteor-trail 
echo is therefore related to an average of the winds in the 1 km vertical 
projection interval. Because of this spatial averaging, a meteor radar is 
not capable of producing the vertical wind structure with as much detail as 
rocket releases, 

A relation between height accuracy and range and elevation angle 
accuracies is found by differentiating equation 2.2: 

j. &h jj-, &h j 

^ = an ^ * H d c 


dh = sine dR - R cose dx. 


(2.3) 
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Then the error in height as a function of errors in range and elevation 
angle is given by 


Ah = sin? AH - F cos? A? (2.4) 

Assuming that A R and A? are normally distributed about zero, the standard 
deviation of h is given by 


= [(sin? a R ) + (R cos? o r ) ] 


1/2 


(2.5) 


where o, , o D , and a are the standard deviations of the height, range and 

rl n l, 

elevation angle measurements . The maximum acceptable value for was 
chosen as .5 km. This value is well under the Fresnel averaging interval, but 
large enough to be realistically achievable. 

Taking R as 140 km and £ as 45°, Table 2,2 gives some representative 
accuracies for range and elevation angle which produce the required height 
accuracy. The first and third table entries concentrate all of the error in 
the range and elevation angle error terms, respectively, of equation (2.2) > 
while the second table entry corresponds to dividing the mean square error 
equally between each of the two error terms. 

2.3.2, Azimuth accuracy requirements . In order to discriminate the 
spatial variation of short-period gravity waves, a horizontal resolution of 
10-100 km is required [ Lindzen ^ 1969]. This is obviously a very weak 
requirement relative to height accuracy; nevertheless a quantitative 
analysis is given below. 

The horizontal displacement from due north, x, is given by 


x = R co ss sinij/ 


( 2 . 6 ) 
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Table 2.2 

Accuracies for the elevation angle and range measurements 
which produce the required height accuracy. 


(deg) 

(km) 

a h (km) 

0 

.7 

.5 

.2 

.5 

.5 

.3 

0 

.5 
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where R is the radial range, £ the elevation angle, and ij/ the azimuth angle 
(Figure 2.7). Differentiating equation (2.6) gives the relation between 
horizontal displacement accuracy, and range, azimuth angle, and elevation 
angle accuracies . Because of the high accuracies required of the range 
and elevation angle measurements by the height tolerance, range and 
elevation angle accuracies should not contribute significantly to horizontal 
displacement error. Under this assumption, the following relation is 
derived from equation (2.6): 

£ | R cos£ cos^j a | (2.7) 

where a ^ and cr^ are the standard deviations of x and ip. Using the same 
values for R and £ as in Section 2.3.1, 


a = 140*. 7 cos'|j a, 
x 

= 100 cosijj 

< 100 a, 

a . > .01 a 
V x 


For a horizontal displacement accuracy of 10 km, 


c. must be about .1 radian. 


or about 6°. In comparison with the elevation angle measurement, azimuth 


angle tolerance is 20 times as great. 

2.3.3 Velocity accuracy requirements. Table 2.1 gives some of the 
characteristics of winds in the meteor region. The mean wind velocity is 
less than about 60 msec and the maximum wave amplitude is about 50 msec"'*'. 
The total horizontal wind should therefore be no larger than about 110-120 
msec ■*. This assumption is supported by the experimental data shown in 
Figure 2.3. A meteor-wind system must, therefore, be capable of measuring 
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horizontal winds as high as 120 msec Meteor-radar systems, however, 
measure radial winds. If the vertical wind component is zero, the radial 
wind velocity is always less than the horizontal velocity. Assuming an 
elevation angle of 45°, radial winds are modified by a factor less than 
cos 45°, or about .7. The capability of measuring radial winds up to 
100 msec" 1 in magnitude should therefore be sufficient. 

Wind velocity accuracy must be a small fraction of atmospheric wave 
amplitudes, as mentioned in Section 2.1, Since wave amplitudes of about 
50 msec -1 will be studied horizontal velocity measurement accuracy of 
about 10 msec" 1 is required. Assuming an elevation angle of 45°, and 
an azimuth angle of 45° also, the radial wind is cos 45°*cos 45°, or only 
.5 of the total wind. Wind accuracy must therefore be .5*10, or 5 msec 1 , 
in order to measure winds with a horizontal accuracy better than 10 msec 
across the entire region of interest. 

In summary, a meteor-radar system must be capable of measuring 
radial velocities up to about ±100 msec" 1 with an accuracy of 5 msec 1 , 

2,4 Processing System Requirements. 

The first constraint on the processing system is that it must keep up 
with the required rate of meteor returns. About 1000 echoes per houT are 
needed; therefore, the system must be able to collect and process the raw 
data from the radar system for at least 1000 returns per hour. 

The output of the radar receiving system consists of actual meteor 
echoes and additive noise. The processing algorithms must operate on this 
noisy output signal, and produce the required measurement accuracies for 
azimuth angle, elevation angle, range, and velocity. Because of noise 
inherent in the receiving system, the signal -to-noise ratio of this output 



25 


signal can be no greater than about 40 dB. The processing algorithms are, 
therefore, constrained to require less than a 40 dB signal-to-noise ratio for 
the input signal in order to produce the measurement accuracies required. 
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3. SYSTEM DESCRIPTION 

A large part of the hardware for the Illinois meteor-wind system is 
taken from the Havana meteor -wind system [Hawkins, 1963] . The Havana system 
was operated for 13 years as a joint project of the National Bureau of 
Standards, Harvard University, Smithsonian Astrophysical Observatory, and 
the National Aeronautics and Space Administration. When the Smithsonian 
program ended in 1971, the system hardware was transferred to the 
University of Illinois. 

The arrangement of the Illinois meteor- wind system differs extensively 
from that of the Havana system. A major emphasis in the Illinois system 
is on moving as much of the control logic as possible out of the system 
hardware and into software, in order to achieve improved reliability. 

This is made possible by the availability of a medium-scale computer, 
capable of controlling the meteor radar; the computer has a good performance 
record as a programmed controller for a partial reflection system and 
a digital ionosonde. 

This chapter describes the function of each software and hardware 
element in the system. In addition, the connection and interaction of 
all the elements is considered. The system hardware elements are first 
characterized; then the interconnection of these elements is given. The 
hardware- software interface and data-frame organization are described 
next. In the last section of this chapter, the software processing 
subsystems are described. 

3 . 1 Hardware Subsystems 

The system hardware is divided into three subsystems for ease of 
description: Transmitting, Receiving, and Processing. In the following 

sections, each is considered in turn. 
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3.1.1 Transmitting. The Illinois meteor-wind system uses a Continental 
Electronics model PO - 830 transmitter, operated at a fixed frequency of 
40.92 MHz. This transmitter is pulse-modulated, and rated for a peak power 
output of 5 Mw. The nominal average power rating is 20 kw, and maximum 
average power 40 kw. The pulse width is constrained by the modulator to 

the range 3 to 100 ysec. In order to maintain the average power level 
at 20 kw, the corresponding pulse repetition frequency (PRF) is limited 
to the range 1330 to 40 Hz. The maximum signal bandwidth allowed by 
the transmitter license is 200 kHz. A pseudo-Gaussian pulse- shape is 
generated by the modulator in order to satisfy this requirement. 

The transmitting antenna for the meteor-radar system is a twelve 
element dipole array, designed and constructed by Aeronomy Laboratory 
personnel. The elements are placed parallel in a single row, spaced by 
.537X and elevated above the ground by . 354X. The main lobe of the 
antenna is oriented due north at an angle of 45° with the horizon. The 
antenna 3-dB beamwidth is 11.7° in elevation and 53° in azimuth at 100 km; 
the gain over a half-wave dipole is 18 dB. The antenna was designed to 
handle a peak power in excess of 5 Mw, 

During both the design and and construction phases, provision was made 
for increasing the gain and horizontal directivity of the antenna by 
adding a second 12-element array beside the first. The antenna will be 
modified in this way if the initial data collected indicates such a change 
to be useful . 

3.1.2 Receiving. In the Illinois System the angles of arrival of 
a radar echo from a meteor trail are measured by phase comparison of the 
received signals at spaced antennas [Lee and Geller , 1973]. The ground 
placement of the receiving antennas is shown in Figure 3.1. Antennas 
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Figure 3.1. Receiving antenna ground placement. 
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A i and 4^ are spaced by one radar wavelength, and form the azimuth pair. 
Antennas A ^ and A ^ are also spaced by one wavelength, and form the 
elevation pair. In addition, antennas A ^ and A^> spaced by 20 wavelengths, 
serve to refine the elevation angle measurement in order to satisfy the high 
elevation angle requirements. A 20-wavelength separation is used because 
of the requirement that elevation angle be measured with accuracy 20 times 
that of azimuth angle. 

Four Telerex Model CY-13-40-42 Yagi antennas are used in the receiving 
system. These 13-element antennas have a gain of 15 dB. The 3-dB 
horizontal beamwidth is 30° , and the vertical beamwidth 64°. The antennas 
are mounted with phase centers at the points 4^ through A 4 on Figure 3.2. 
They are oriented due north, and elevated at 45°. 

Four ARI PR 40-A receivers are used in the system. Each receiver 
has a fixed gain of 100 dB, a bandwidth of 180 kHz, and a noise figure 
less than 4 dB, and includes an envelope detector and two phase-comparator 
output channels. The envelope detector operates in two modes: linear 
and logarithmic. In the linear mode, it has a dynamic range of 30 dB, 
and the signal output at the top of the linear region is 1 V. The phase- 
comparator output channels produce signals proportional to the sine and 
cosine of the phase difference between the received signal and a reference 
input. The dynamic range and output of the phase- comparator channels 
are similar to those of the envelope detector. 

Each receiver is connected to an antenna through a length of coaxial 

cable. The signal loss through the cable is 6 dB. For a sky-noise 
4 

temperature of 10 K, the rms noise at the output of the receiver is 

.2 V. The receiver output noise with the input connected to a 50-J2 resistor 



is .05 V. Then about 6 dB of attenuation is required at the input to each 
receiver, in addition to the 6 dB cable loss, in order to scale the 
input signal and make the maximum use of the receiver dynamic range. 

3.1.3 Processing. The meteor-radar system is supported by a 
Digital Equipment Corporation Model PDP-15 digital computer system. The 
central processor of this system is an 18-bit machine, with 32 k words 
of .8-psec core. An extended arithmetic element provides hardware multiply 
and divide. The machine has a real-time clock, and a data channel for 
high speed input /output . 

Several basic peripherals are part of the system. Four DECtapes 
provide medium-speed bulk storage, and two disks are used for high-speed 
on-line storage. A 50-character -per-second paper- tape punch and a 300- 
character-per-second paper-tape reader are used for low-speed storage. 
Inputs are from a Teletype, Model KSR 35, which prints 10 characters per 
second, and a Hazeltine Model 2000 cathode ray terminal, which operates at 
870 characters per second. 

Analog input is provided by a Hewlett-Packard Model 5610 10-bit 
A/D converter. The device performs 10 5 conversions per second, and 
has a sampling aperture of .05 ysec. The A/D converter has an accuracy 
of *1/2 least significant bit, and uses two ' s-complement output coding. 

A 16-channel input multiplexer selects a signal line to be sampled. Input 
signals in the range from -1 V to .998 V are mapped onto integers in the 
range from -512 to 511. The digital output of the converter enters the 
system through a direct-memory-access interface. 

Digital input and output is also available. An interface allows 18- 
bit words to be sent or received by the processor through the accumulator. 
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The signals for the device are TTL logic levels. The maximum transfer 

4 

rate for input or output is about 10 words/sec. Access to the input 
and output signals is provided by two connectors mounted on the front panel 
of the computer. 

A Background/ Foreground software monitor manages the system. The 
monitor supports assembly programming, and FORTRAN IV compilation. All 
system peripherals are supported by the monitor, and all system programs 
are disk-resident for fast access. 

An Automatic Priority Interrupt (API) system controls the asynchronous 
operation of the various peripheral devices and software tasks. Of the 8 
priority levels available, the lowest four are assigned to software systems 
and the highest four to hardware devices. The software priority levels are 
used to set up a priority queue for processing real-time data, without 
inhibiting hardware interrupts from peripheral devices. 

3.2 Hardware Interconnection 

There are 12 receiver channels: two phase comparators and one envelope 
detector output for each of four receivers. However, there are only 9 
unique information channels, since the signals in the four envelope-detector 
channels are essentially the same. The envelope-detector channel from 
the receiver connected to antenna 1 is used, and the other 3 are ignored. 

The information from the 8 phase-detector and one envelope*»detector channels 
is sent to the A/D converter through the multiplexer. 

In order to simplify the hardware system as much as possible, syn- 
chronous control is used. A master oscillator generates the pulse- 
repetition frequency, and the transmitter, multiplexer, and A/D converter 
are all synched to this oscillator (Figure 3.2). The operating sequence 
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is as follows: A pulse from the master oscillator causes an R F pulse 

to be transmitted, and also causes one channel to be selected by the 
multiplexer for input to the A/D converter. After a delay, the A/D 
converter begins sampling the selected channel, and sending the digital 
values corresponding to the input voltage at the sampling instant to 
the computer. The delay period is used to ignore echoes from objects 
too close to be meteor trails. After the required number of samples 
is taken corresponding to the radial region of interest, sampling ends* 

The sequence begins again with the next pulse generated by the master 
oscillator . 

There are several consequences of this scheme. Note that the 
sequence is not affected by the data in the receiver channels. Whether 
or not an echo is present, the sequence continues. Secondly, for each 
pulse transmitted, only one channel is sampled. During the time when an 
echo is present, only 1/9 of the useful data is collected for each pulse 
transmitted. A third consequence is that an upper limit is placed on the 
useable PRF range. A pulse cannot be transmitted until sampling by the 
A/D converter is ended. 

The delay period and sampling interval are selected by considering the 
region of interest (Figure 3.3). All meteor echoes come from the 80-120 km 
height interval. Considering the transmitting and receiving antenna ele- 
vation beamwidths, the range of elevation angles for received echoes is 
about 45° ±15° . Then the minimum time after pulse transmission for which 

a meteor-trail echo can occur, T . , is 

> min’ 

_ 2 * 80 
min ~ sin 6CT • a 


= 616 ysec 



RFAM AVIS 



MITTER 


Figure 3.3. Meteor echo region. 
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where a is the speed of light. The factor of 2 in the numerator comes 

from the round-trip path of the radar pulse. The maximum time of interest 

after pulse transmission, T , is 
r ’ max 

2-120 
tV = . .HUM— I . 

max sin 30 • c 

= 1600 nsec 

The maximum usable PRF is simply T * , or 625 Hz. The A/D converter 

max 

sampling period is the difference between T and T . , or 984 psec. 

ilidA JUXIX 

This time period, however, must be a multiple of the A/D converter 
sampling interval, which is 10 psec. Times and T max are, therefore, 

chosen as 610 and 1600 psec, in order to cover the region of interest 
and still satisfy this requirement. 

The various hardware elements must be interfaced. The reference 
signal into each of the receivers is simply the RF signal used by the 
transmitter. The receiver-output channels are connected directly to 
the multiplexer and the A/D converter, since the linear region for the 
receiver channels corresponds to the non-saturated region for the A/D 
converter (i.e., ±1 V) . A TTL- level signal is needed to sequence multi- 
plexer channels and also to initiate sampling by the A/D converter. The 
transmitter, however, requires a 5-V input pulse for triggering. The 
PRF generator delivers 5-V pulses. Therefore, an interface is needed 
between the PRF generator and the multiplexer. The delay element must 
accept a 5-V pulse as input, wait for 610 usee, and then deliver a TTL- 
level signal to the A/D converter. This signal is the blanking pulse; 
it must have a "0” active state, and must remain in the "0" state for the 
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duration of the sampling period. A timing diagram of the required signals 
is given in Figure 3.4. 

3 . 3 Data Frame Organization 

In the data- col lection system, a hierarchical data structure is 

used, consisting of elements, lines, and frames. A line is the sequence 

of digital values sent to the computer by the A/D converter during each 

sampling period. The individual values that form a line are elements. 

The sampling period covered by one line is 2 ’ max - J ' min , 990 usee. Elements 

are taken at 10-psec intervals, starting at T . and ending at T A 

mm & max * 

line is, therefore, a sequence of 100 elements. A frame consists of the 
number of lines required to measure all the parameters of one meteor-trail 
echo. Since there are 9 different receiver channels, at least 9 lines 
are needed to form one frame; in fact, more are required. 

In order to determine the number of lines needed to make a frame, 
the requirements of the various measurement algorithms are considered. 

The echo-detection and range algorithm together require one line. Echoes 
are detected by searching a line taken from the receiver envelope- detector 
channel for elements with magnitude above a noise threshold; the range 
algorithm operates by considering the position in the line of these ele- 
ments. The phase-comparison algorithm for antenna pairs requires four 
lines. These lines are taken sequentially from the sine and cosine channel 
first one receiver, and then the other. In addition the Doppler frequency 
is needed. The Doppler frequency algorithm measures the change with time 
of the phase of the received signal at any of the four antennas. It re- 
quires that a total of four lines be collected alternately from the sine 
and cosine channels of any one receiver. Detailed descriptions of these 
algorithms are given in the next chapter. 
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Figure 3.4. Delay element timing diagram. 



Each frame must contain the data needed by a number of algorithms. 

Echo detection and range measurements must be made twice per frame to 
give good noise rejection. Three phase comparisons are needed: one each 

for the signals from the azimuth, coarse elevation, and fine elevation 
antenna pairs. The Doppler shift of each echo must be measured twice 
per frame, since the phase-comparison algorithm uses the Doppler measurement. 
The echo-detection algorithm uses one line, the Doppler and phase-comparison 
algorithms use four lines each. Then apparently 2 + 12 + 8, or 22 lines 
are needed per frame. By carefully arranging the order of the measurements 
in the frame so that some of the data lines are used by more than one 
algorithm, this number can be decreased. 

The order of measurements is chosen with certain guidelines in mind. 

The two echo detection and range measurements should be well separated in 
the frame so that any noise present in the first line will be decorre- 
lated from that in the second. Also, the fine-elevation phase comparison 
should be preceded and followed by a Doppler measurement, since this 
phase comparison is the most critical, and the phase-comparison algori thm 
accuracy is strongly related to Doppler accuracy. 

Using S. and C. to symbolize lines taken from the sine and cosine 

i / 

channels respectively of the receiver connected to antenna i, and V as 
the video output of the envelope-detector channel. Figure 3.5 shows a 
sequence of 16 lines which hold the data for the required measurements 
in the proper sequence. The total data frame is, therefore, 16 lines of 
100 elements, for a total of 1600 elements; one data frame is collected 
for every 16 pulses transmitted. Since exactly 16 lines are needed, 
the 16-channel multiplexer can be used to good advantage. Each receiver 
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channel is connected to the multiplexer channel corresponding to its 
position in the 16-line sequence shown in Figure 3,5. Note that some 
of the receiver channels are connected to more than one multiplexer 
channel. With the receiver channels connected to the multiplexer in this 
way, the multiplexer simply sequences through each of its input channels 
in order. One pass through all the input channels produces one data 
frame. 

3.4 Data Processing Subsystems 

The data processing system (Figure 3.6) is made up of three parts: 
collection, search, and reduction. The collection program controls the A/D 
converter and runs at the highest priority level. This program operates 
synchronously with the hardware system. Lines are collected from the A/D 
converter until a full frame is available; then, control is passed to 
the echo-search program at a lower priority level. Collection continues 
even while the echo-search program is running, since the priority level 
of the collection program allows all other programs to be interrupted 
whenever necessary. The echo-search program scans the two video-output 
lines of the frame. If an element above the noise threshold is found, 
a reduced frame is built; otherwise the frame is discarded. A reduced 
frame is made up of certain essential elements taken from the total 
frame, which characterize the echo. Each reduced frame is pushed onto 
a first-in, first-out stack for temporary storage. The reduction program 
takes a reduced frame as input data, and produces the polar coordinates 
and velocity of meteor- trail echoes. 

There are three main data storage areas. The first is the frame 
storage buffer. Data are written into this region by the collection 
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Figure 3.6. Data processing system interconnection. 
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program, and read by the range search and echo detection subsystem. 

Since concurrent collection and echo search is needed, a double buffer 
is used. That is, collection and search are done alternately in two 
buffers. Since one frame has 1600 elements, a 3200 word storage area 
is needed. The second major storage area is the first-in, first-out stack. 
The echo detection and range program produces a 22-element reduced frame, 
which is stored temporarily in this stack, until the data reduction 
program is ready to accept this information as input data. The stack 
has 220 words, enough memory to store 10 reduced frames. The third 
major storage area is the system status table. This table is only 10 
words long, but it contains all the information needed for communications 
between the various subsystems. The function of each word is given in 
Table 3.1. 

These three storage blocks are Fortran Named Common storage areas; 
the frame buffer is called FRAME, the stack is named STACK, and the 
status table is named STATUS. Only STATUS is initialized by a BLOCK 
DATA subprogram. 

Program modularity is an important consideration in real-time 
data processing systems, since changes are often made. By breaking the 
system into small modules, each performing one particular function, 
modifications are easily made, and unexpected side effects from small 
changes are avoided. With this in mind, two rules are imposed on 
the software coding: 

1) No single program or subprogram shall have more than 50 


executable lines. 



Table 3.1 




Word definitions for the system status table 


Word 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Function 

Time in seconds since start of collection 
Pointer to current line being filled 
Pointer to current frame being frame being filled 
Pointer to last reduced frame written on stack 
Pointer to last reduced frame read from stack 

y not yet assigned 
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2) No program or subprogram shall have more than one 
entry point . 

3.4.1 Data collection. The data collection program is an interrupt 
service routine for the A/D converter. Each time a line is completed, the 
processor is interrupted and the collection program is executed. After 
the collection program is finished, the interrupted program is restarted. 
The collection program performs two operations. First, the A/D converter 
is set up to collect the next line. Then, if a frame has just been 
completed, the echo search and range program is started. Figure 3.7 

is a flow chart of the collection program. In this figure, J is a 
pointer to the current line. 

3.4.2 Meteor echo search . The maximum-amplitude element in a video 
line, along with the elements immediately preceding and following, is 
called the range set. If the location in a video line of the maximum 
amplitude element is position £, then the range set for this line is 

the group of elements i-1, i, and £+1. The meteor echo search begins 
with selection of the range set in the first video line. Next, elements 
i-1, i, and i+1 are selected from the second video line. A threshold 
test is done next. If the element in position i of both video lines 
is above the noise threshold, and neither is saturated (that is, neither 
has a value of 511 or -512) then a usable meteor trail echo is present. 

If no such echo is present, the frame is discarded. If a trail echo 
is present, however, the following actions are taken. The six elements 
taken from the two video lines are pushed onto the stack, along with i, 
the position of the center sample in the first video channel range set. 
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Figure 3.7. Data collection routine. 
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Each of the elements in position i of the 14 sine and cosine lines are 
also added to the stack. Finally, the current time in seconds relative 
to the beginning of data collection is put on the stack. Altogether, 

22 samples are recorded. These 22 samples form one reduced frame. A 
flow chart of this algorithm is shown in Figure 3.8. 

3.4.3 Data reduction. The data-reduction subsystem takes reduced 
frames from the stack, and calculates the polar coordinates and radial 
velocity of meteor-trail echoes. The program uses all the data reduction 
algorithms which are described in the next chapter; each algorithm is 
a subroutine. 

The reduction program itself does four things: it interfaces to 

the stack, formats the data for the algorithm subroutines, performs 
consistency checks on the redundant frame measurements, and writes the 
results of the reduction onto magnetic tape. 

The range algorithm is applied to the data from both video lines. 

A consistency check is made on the results; if the results do not agree 
within the required range accuracy of 700 m, the reduced frame is rejected 
and no further processing is done. If the consistency test is passed, the 

Doppler algorithm is applied to both Doppler data sets. Again, a consistency 
check is made; the results must agree to within 5 msec . The various 
phase comparisons for the antenna pairs are then made, after which the 
phase differences for the antenna pairs are combined with antenna spacing 
information to give the angles of arrival. Finally, all of the results 
are written onto magnetic tape. A flow chart of the reduction subsystem 
is shown in Figure 3.9. 




Figure 3,8. Meteor echo search routine. 





Figure 3.9. Data reduction routine. 
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4. MEASUREMENT ALGORITHMS 

This chapter describes in detail each of the four measurement 
algorithms which are used in the processing system. The measurement 
accuracy and magnitude limitations of the algorithms are also con- 
sidered. The effect of the PRF and pulse width on measurement accuracy 
is examined, and the pulse width and PRF which optimize measurement 
accuracy are found. 

4.1 Range Determination by Curve Fitting 

In meteor radar, range is deduced from the time difference between 
the instant of pulse transmission and the instant of echo reception. 

The instant of pulse transmission is that instant during pulse transmission 
when the maximum voltage is present at the transmitting antenna terminals, 
while the instant of echo reception is that instant during echo reception 
when the maximum voltage is present at the envelope detector video output. 
Range is calculated by multiplying this time difference by e/2. Of 
course, any system delays introduced by the receiver or interconnecting 
cables must be subtracted from the time difference before the range 
calculation. 

In the Illinois system, the time difference is measured in three 
parts. First, there is the 610-ysec duration of the blanking pulse 
(Section 3.2) established to reject radar clutter. Second, there is 
the time between collection of the first element in the video line, at 
the end of the blanking pulse, and the instant of sampling of the largest 
element in the line. The elements are samples taken at equal 10-ysec 
intervals, so if the ll/th element is the largest, this second time 
interval is 10(iV-l)usee. The instant the largest sample is taken does 
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not correspond exactly with the instant the largest voltage is present 
at the receiver video output; so a third part of the time difference, 
x (Figure 4.1) corresponds to the interval between the largest voltage, 
and the instant the largest element is sampled. Range is calculated 
according to the expression 

a = (610 + 10 (ff - 1) + x) * 1.5 (4.1) 

where R is the range in kilometers, and x is in usee. For a symmetric 
pulse, t is in the range ±5 usee. If x is assumed zero, the range may 
be in error by ±750 m. Recall from Section 2.3.1 that the required 
range accuracy is ±500 m. Therefore, some method is needed for measuring 
x. This method is the range algorithm, which uses a process of curve 
fitting. Two aspects of the sampling of the video output must be 
considered before fitting a curve to the sampled elements. These aspects 
are the finite aperture time and sampling rate of the A/D converter, 
and are discussed in the following paragraphs. 

The finite aperture time of the A/D converter causes an integration 
of the input signal; however, the video-output voltage is slowly varying 
relative to the .05-usec aperture time of the A/D converter since the 
receiver has a 180-kHz bandwidth. Then the output of the A/D converter 
is a good approximation to the input voltage at the instant of sampling, 
and the finite aperture time has little effect on sample accuracy. 

The radar system has a characteristic pulse shape, defined as the 
time record of the voltage at the receiver video output during reception 
of an echo from a hard reflector. This shape is determined only by the 



At = 10 iLsec 


Figure 4.1. Typical receiver output pulse, with sampling superimposed 
The quantity r is the offset between point of maximum on 
the pulse, and highest sample, and A A and A are the 
values of samples occurring during pulse reception. 
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transmitted pulse shape and the receiver bandpass, and is the convolution 
of the receiver impulse response with the transmitted pulse. Of course, 
convolution in the time domain is simply multiplication in the frequency 
domain. Both the receiver bandpass and the transmitted pulse are nearly 
Gaussian. Using the following Fourier -transform pair, the characteristic 
pulse can be calculated from the receiver bandpass and the transmitted 
pulse : . r \ 


2tt) = (-) 


1/2 


exp (-st )«-*■£(/) = exp 


■i£ 


( 4 . 2 ) 


The transmitted pulse is transformed to the frequency domain and then 
multiplied with the receiver bandpass. The product is then transformed 
to the time domain, giving the characteristic pulse. For a Gaussian 
receiver bandpass and transmitted pulse, the characteristic pulse is 
also Gaussian. By comparing the characteristic pulse with the transmitted 
pulse, the shape modification due to the receiver is found. The ratio 
r of the characteristic pulse length to the transmitted pulse length, 
and the ratio r of the characteristic pulse maximum amplitude to the 
transmitted pulse maximum amplitude, are given by 


« [1 + (ir BL )~ 2 ] 
= [1 + (t rSI)' 2 ] 


1/2 

- 1/2 


( 4 . 3 ) 

( 4 . 4 ) 


where B is the receiver bandwidth and L the transmitted pulse width. 
The level used to define the bandwidth and pulse width is arbitrary; 
however, the same level must be used in both definitions. Also, 
consistent units are required, i.e,: Hz and sec, MHz and psec. 
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The sampling rate of the A/D converter is 10^ sec~\ so if the 
characteristic pulse is less than 10-psec wide, sampling may miss the 
pulse entirely. For a pulse width of 10-20 psec, at least one sample 
will always occur during pulse reception; for a width of 20-30 psec, 
at least two will occur; for a pulse width greater than 30 psec, three 
or more samples will occur during pulse reception. There is no way 
to compute t if the pulse width is under 20 psec. 

A technique for finding t is available if three criteria are met: 

1) The characteristic pulse shape is known. 2) The characteristic pulse 
width is such that more than one sample is taken in its duration, 3) The 
received echo is from a hard reflector. Under these conditions, the time 
record of the video output can be reconstructed from the samples taken 
during pulse reception. This reconstruction amounts to fitting the 
known characteristic pulse to the samples. Once the time record of the 
video output has been reconstructed, the quantity t is readily available. 
Since samples need only be above the noise level to be used in the pulse 
fitting, the width of the characteristic pulse is defined as the time 
between 20-dB points. The pulse fitting method used for pulse 
widths in the range from 20-30 psec takes the ratio of the maximum 
element to the higher of its two neighbors. The method used for pulse 
widths greater than 30 psec fits a parabola to the maximum element and its 
two immediate neighbors. The meteor-wind system uses the parabola method, 
for reasons given in Section 6.1. Both methods are described below, 
however, for completeness. 
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4.1.1 Parabola method. The parabola method assumes that the 

shape of the characteristic pulse between 20-dB points can be modeled 

successfully by a parabola. For the Gaussian- shaped pulse that the 

meteor -radar transmitter produces, this assumption is reasonable. The 

procedure for finding t is to first find the coefficients of the 

2 

function f(x) -ax + bx + a, which passes through the maximum element, 
and each of its immediate neighbors. The first derivative f’(x) is 
found; t is equal to the value of x for which /' (r) = 0. Elements A^, 

A 2 > and A ^ are the three points to which the parabola is fitted (Figure 
4.1). Relative to the center element, A ^ and A ^ are taken at x = -10 
and x - 10, respectively. Then three equations can be written. 


a(-10) 2 

+ 

KF 

1 

h— 1 

o 

+ o = A 

a(0) 2 

+ M0) 

+ a = A 

a(10) 2 

+ b (10) 

+ a ~ A 


These equations are then simplified and solved. 
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(4.5) 


The first derivative is calculated and set to zero: 
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f' (a;) = 2ax + b 

r (t) = o 

_ -& 

T " 2 a 

Quantities a and b are substituted from equation (4.5) : 

, V A 3 

' W 2A 2 


(4.6) 


(4.7) 


Then the offset x, in ysec, is found according to equation (4.7) . 

Note that the parabola method does not require prior knowledge of the 
characteristic pulse shape. Even if the pulse shape is distorted by 
reflector irregularities, the method still gives a reasonable approximation 
to t. The only limitation is the modeling of the pulse shape between 
20-dB points by a parabola. For a Gaussian-shape characteristic pulse, 

30 usee wide, the maximum error in x due to parabola modeling is 1 ysec. 

For no-noise conditions, then, the parabola- fitting technique gives a 
range accuracy of ±150 m; this accuracy degrades gracefully with pulse 
distortion and atmospheric noise. 

4,1,2 Ratio method . If the characteristic pulse is 20-30 ysec wide, 
only two samples will be taken during pulse reception; finding x by a 
parabola fit is not possible. If the shape of the pulse is known, 
however, the ratio of the two samples can be used to imply x. For 
any given pulse shape, a table can be formed which relates x to the 
ratio of the two samples. For a symmetric pulse, the points of pulse 
maximum must occur somewhere in the closed interval between the maximum 
sample, and the midpoint between the two samples; that is. Os |x| s 5 ysec. 
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Table 4.1 is developed by assuming that a Gaussian pulse is transmitted, 
and that the receiver has a Gaussian bandpass. For various receiver 
bandwidths, the ratio for several values of t is given. Negative values 
of t are not given, because in this case the table entry is just the 
reciprocal of that for positive t. 

The accuracy limitations of this technique are twofold. First, 
the largest practical size for the table is about 6 elements as 
shown, implying an accuracy of 1 ysec in t. Second, pulse distortion 
and noise introduce large errors. For these reasons, the ratio technique 
is limited to an accuracy of about ±2 ysec. The corresponding range 
accuracy is ±300 m. 

4 . 1,3 Conclusion. For the range algorithm, the parabola and ratio 
methods both give satisfactory results. The parabola method, however, is 
less susceptible to errors from pulse distortion and low signal-to-noise 
ratio then the ratio method and seems to be the optimal technique. For 
a given transmitter average power, a wide pulse implies a low PRF. It 
is shown in the following subsections, however, that the Doppler and 
phase- comparison algorithms are most accurate when a high PRF is used. 

The selection of the optimal pulse width and PRF must be deferred until 
the accuracy of the entire system has been considered. The final 
selection is made in Section 6. 

4.2 Doppler Measurement 

Each system receiver has a phase comparator. Phase comparisons 
axe made between receiver IF signals and a 5.5 MHz signal derived from 
and phase- synchronous with the transmitter oscillator. The frequency 
of the IF signal in each receiver is 5.5 MHz plus the Doppler shift of the 
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Table 4.1 

Ratio of the highest sample to the higher of its two neighbors for 
offset x in the range 0-5 usee, and receiver bandwidth in the range 
20-150 kHz. Transmitted pulse width is 5 usee. Transmitted pulse 
and receiver bandpass are Gaussian. 


Offset (microseconds) 




0 

1 

2 

3 

4 

5 


20 

1.43 

1.33 

1.24 

1.16 

1.08 

1.00 


50 

4.60 

3.39 

2.50 

1.84 

1.36 

1.00 

Bandwidth 


75 

10.23 

6.42 

4.04 

2.53 

1.59 

1.00 

(kHz) 


100 

17.22 

9.75 

5.52 

3.12 

1.77 

1.00 


150 

29.65 

15.05 

7.64 

3.88 

1.97 

1.00 
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received signal. The phase comparators produce two outputs, signals 
proportional to the sine and cosine of the phase difference between the 
IF signals and the reference input. The time derivative of the phase 
at any of the receivers is the Doppler-frequency shift of the received 
signal. A method for calculating this time derivative from elements 
taken from the sine and cosine outputs of one of the phase comparators 
during echo reception is needed. This method is the Doppler algorithm. 

The outputs of any phase comparator at some instant t can be written 

A cos (6 + ait) 

A sin (0 + wt) 


where A = amplitude 

a) - Doppler frequency 
0 = phase of signal at antenna 

Doppler frequency is measured as follows: at time t measure the 

cosine-output signal; at time measure the sine-output signal; 
at time ty measure the cosine-output signal again. These measurements 
are labeled M y My the Doppler frequency is given by the relation 
(Appendix I) 


u 


sin^Z 1 

At 


( 4 . 8 ) 


where 


Z* = 


M -M 

1 3 

2 


At = = £3 “ *2 = P u ^ se repetition period 
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= A sin(0 + (otp 

M 2 - A cos(0 + 

= A sin(0 + wtp 


Quantities AL, M 2 , and A/, are input through the A/D converter, 
causing some quantization error in each. The worst- case error in Z* 
due to quantization of Afj , and Af is calculated in Appendix II, 
with the result 


| AS ' | < 


1/2 + sinjojAt 


Wo 


(4.9) 


At this point, it is covenient to define a quantity a, the Doppler 
angle, as the phase change due to Doppler shift at any receiver in each 
pulse -repetition period. That is, 


ct = u)A t 

Then the expression for AZ ’ is rewritten 


AZ' 


1/2 


+ sin 




kL 


(4.10) 


Note that the denominator in this expression can be small; large 
errors in Z' result. If a fourth measurement, M ^ = cos (6 + ^) , is 

taken, an alternative method of finding sin toAt is available. A 
trigonometric argument similar to that in Appendix I gives 


W 4 -W 2 

2M- 


Z" = 


= sin toAf 


(4.11) 
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Quantization error in this expression is given by 


AZ" 


1/2 + sin a 




(4.12) 


Both AZ ' and AZ" are of the same form, and differ only in the denominator. 
The method for finding sin coAt which gives the least error is, therefore, 
the method which has the largest- denominator. Two new quantities, Z and A/, 
are defined: 

if \M 2 \ > \M Z \ , M = M 2 and Z=Z* 
if |M 3 1 £ \M 2 \ , M = M and Z=Z" 

The Doppler algorithm consists of calculating Z and finding a, where 


a = sin Z 


(4.13) 


Then the expression for the error in Z is 


| AZ | < 


1/2 + sin | a 


M 


(4.14) 


The question arises: what is the smallest value for the magnitude 

of M as a function of ot? This question is considered in Appendix III, 
with the result 


1*1 > 


A sin(Tr/4 - |^j) 


Then equation (4.14) is rewritten 


AZ | < 


1/2 + sin a 


A sin(ir/4 - ||-|) 


(4.15) 
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The velocity of a reflector relative to the transmitter is related to the 
Doppler shift by 


wX _ a\ 

4 it 4ir At 


(4.16) 


Then errors in velocity are related to errors in Doppler angle by 


Ay 



(4.17) 


The quantity Aa is gotten by differentiating equation (4.13) 


But Z = sina, so 


da = 


d Z 


A-z* 


Aa = 


_ AZ 


cosa 


Aa | < 


1/2 sin |a j 
A sin(Tr/4 It!) cosa 


< /(a) 


(4.18) 


(4.19) 


(4.20) 



(4.21) 


Inspection of equation (4.20) reveals that /“(a) is a monotonically increasing 
function; a 1 is monotonically decreasing. Then /(a) a"" 1 has a minimum 
(Figure 4.2). Numerical evaluation of this function shows that this minimum 
occurs for a = .446 radians. 
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From the definition a = coAt, the optimal technique for measuring 
oj seems to be to vary At, the time between transmitted pulses, until 
ot = .446. Then u * .446/At. This method is impractical because the 
lifetime of a meteor trail is finite. For low velocities. At may become 
larger than the trail lifetime. A practical approach is to choose At 
such that a = ,446 for the maximum velocity of interest. Since this 
maximum velocity is about 100 msec -1 , 

(*) A 
. max 

v = — j— 

max 4 it 


100 * 4TT n, . j -1 
w = — ~ a — = 171.4 rad. sec 
max 7 . 33 


Then 


a - u At 
max 


446 

At = — 4 ~ * 2 * 602 msec 


The optimum pulse-repetition frequency for a single measurement of 
Doppler frequency is then 384 Hz. If more than one measurement is to 
be taken, another factor must be considered. For any given reflector 
lifetime, the number of measurements that can be taken is proportional 
to the PRF. The improvement in accuracy of a measurement due to data 
averaging is proportional to the square root of the number of values 
averaged. Equation (4.21) is modified to 


Av 


t ~ 


££al . _H 

a v^PRF 


C4.22) 
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But a is inversely related to the PRF, so that 

to . ~m. va in 

a 

Ay' ~ /(a) a 1 y (4.23) 

- 1/2 

When the quantity f[a)a is plotted, the minimum occurs for a = .227 rad 

(Figure 4.2). Following the procedure used for the single measurement 

case, the optimum PRF is 755 Hz. It is important to note that the function 
- 1/2 

f(a)a is fairly flat near the minimum. The implication is that the 
advantage gained by multiple measurement is very nearly cancelled by the 
loss in accuracy due to a higher PRF. One concludes that any PRF in the 
range 400-1700 Hz, corresponding to a in the range .1-.4 rad, is acceptable. 

The single-measurement accuracy is calculated using the following 
parameters : 

y = 100 msec"'*' 

PRF = 400 Hz 

47ry , , -1 

to = — — ~ 171.4 rad. sec 

a = aiAt = .4285 rad 



100 


435 
' A 

Since A, the signal amplitude, can range from 0-511, a reasonable value 
is 100. Then | Ay ] <4.85 msec - * - . 



4 . 3 Phase Comparison 

In order to calculate angles of arrival, phase comparisons between 
pairs of separate antenna-receiver sets must be made. Such phase 
comparisons are required for measurement of coarse elevation, fine 
elevation, and azimuth angles. The phase-comparison problem can be 
easily understood in the light of the following analogy: given two 

clocks which cannot be viewed simultaneously; how can their indicated 
times be compared? If one observes the first, and then, a moment 
later, observes the second, both clocks have moved forward by some 
time increment. Obviously, one must use a stopwatch in order to 
know the displacement which occurred between observations. In the 
case of meteor radar, clocks one and two correspond to phase at two 
receivers. The phases cannot be observed simultaneously because 
receiver outputs are multiplexed sequentially (Section 3.3). The 
required "stopwatch" is simply the Doppler frequency; the phase 
increment between measurements is the Doppler frequency multiplied 
by the time increment between measurements. This quantity is a, 
the Doppler angle, as defined in Section 4.2. 

The phase comparison requires the following sequential measurements 

1. Sine output of Receiver #1 

2. Cosine output of Receiver #1 

3. Sine output of Receiver #2 

4. Cosine output of Receiver #2 

From the above analogy it appears that there should be some way 
to use these four quantities, together with Doppler frequency and time 
between measurements, for a calculation of phase difference. Such a 
calculation is outlined in Appendix IV. 
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Singularities in this calculation can only occur when all four measure- 
ments return the value zero. This case is possible only when 

wA t = tt/2 

w = Doppler frequency 
At = pulse-repetition period 

Since both the sine and cosine of the phase difference are given by 
the calculation, phase differences can always be calculated within ambi- 
guities of multiples of 2tt radians. This ambiguity is eliminated by 
antenna placement. 

Accuracy of the phase -difference measurements is, to a first order, 
the same as the accuracy of measurements of the Doppler angle as defined 
in Section 4.2. For this reason, phase-difference measurement is accurate 
to at least within .03 radians. 

4.4 Angles of Arrival 

In the Illinois meteor-wind system, the angles of arrival of an 
echo are obtained using four spaced receiving antennas in an inter- 
ferometric mode. To determine the angles of arrival, the phase 
differences between pairs of antennas are measured. The symbols 

6 , . and d . . are used to represent the phase difference and ground 

I'd I'd 

distance between antennas A. and A.. Then the phase differences of 

i 3 

interest are written (see Figures 3.1 and 3.5) 

2tt 

$ 12 M — i — “ cose cosif* , (4.24) 
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2n d 


41 


P 41 


COS? COSlj) 


(4.25) 


2ir ^ 

<j>l 3 = j— cost; sin<|> 


(4 . 26) 


where ? and tl> are the elevation and azimuth angles of arrival , and 
X the radar wavelength. Phase differences and are taken from 
antennas aligned along the beam axis. Distances d ^ and d ^ are 
IX and 20X, respectively. These two phase differences are used to 
find the elevation phase difference, with the widely spaced antenna 
pair acting as a vernier to refine the phase difference determined by 
the closely spaced pair. The azimuth phase difference is simply <t>^, 
which is taken from the antenna pair spaced one wavelength, and oriented 
perpendicular to the beam axis. 

Two quantities, P ^ and P , are defined as 


P 41 


el 2tt d 


41 


31 


az 


2tt d 


31 


(4.27) 

(4.28) 


where 4,. is the refined version of 4,,. Then P 7 and P are the 
T 41 T 41 el az 

normalized phase differences for the elevation and azimuth antenna 
pairs ; 

P^£ = cos? cos ip 

P = cos? sini|> 
az 


(4.30) 
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The angles of arrival algorithm is a three-step procedure. First, the 
wide-spaced elevation phase data are used to refine the close-spaced 
elevation phase information. Next, the elevation and azimuth phase 
differences are normalized. Finally, the normalized phase differences 
are combined to yield the azimuth and elevation angles of arrival. 
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5. DATA REDUCTION SUBSYSTEM CODING AND SIMULATION TESTING 

The data-reduction subsystem contains all the measurement algorithms, 
and operates on the reduced frames produced by the echo- search subsystem. 

In order to verify the proper operation of the measurement algorithms, a 
simulation of the hardware and software subsystems up to data reduction is 
done. That is, a simulation program is used which, given as input data 
the parameters that describe a meteor-trail echo, produces a reduced 
frame. These descriptive parameters include the coordinates and velocity 
of the specular point, and also the amplitude of the echo signal and 
noise at the receiver output. The reduced frames produced by simulation 
are input to the reduction system, and the coordinates and velocity output 
by the reduction system are compared to the echo parameters initially 
given to the simulation program. The simulation is used to test the data- 
reduction subsystem on a comprehensive set of parameters which include 
all of the expected characteristics of meteor-trail echoes. In this way, 
algorithm validity is demonstrated, and the minimum signal-to-noise ratio 
which gives the required measurement accuracies is determined. 

Data-reduction subsystem testing is done in three steps. First, the 
measurement algorithms are used to design the data-reduction program. 

Next, the program for simulating the meteor-radar system up to the data- 
reduction subsystem is written. Finally, the actual testing is done. 

Each step is described in detail in the following paragraphs. 

5.1 Subsystem Coding 

The data-reduction subsystem consists of 6 subroutines, 4 Common Blocks, 
and 2 Block Data programs . Each subroutine performs a single operation 
on the data stream. Three of the Common Blocks are used for communications 
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between subroutines; the fourth block contains data related to the spacing 
between antennas. The six subroutines which perform the data 
reduction are described here; listings are given in Appendix V. 

SDRIV1 is the interface between the reduced frames input to the 
reduction system, and the subroutines containing measurement algorithms. 

It organizes the data for use by the measurement subroutines. Each of 
the other subroutines is subject to certain error conditions, which can 
occur with a noisy echo or a very high Doppler shift. These error 
conditions are sensed within the measurement subroutines, and communicated 
to SDRIV1 through one of the common blocks. 

SDRIV1 groups the input data and calls the various measurement 
subroutines in the proper order. If no error is sensed in any of the 
measurement subroutines, SDRIV1 outputs the measured quantities; otherwise, 
an appropriate error message is printed indicating the nature of the error. 

DOPPLER is the implementation of the Doppler algorithm. Elements 
taken from four of the phase lines are its arguments, and sina is its 
output. One error condition exists. If both the second and third 
arguments are zero, no calculation can be done. That is, if both M 2 
and M ^ are zero from Section 4.2, an error condition occurs. This error 
is flagged as error 21. 

SINCOS takes sina, the output of DOPPLER, and generates the sine 
and cosine of multiples of a as well as a itself. The sine and cosine 
of multiples of a are needed by the phase-difference algorithm. The 
quantity a is needed for determination of the radial velocity. An error 
occurs if the input quantity sina has magnitude greater than unity. This 
condition is flagged as error 41. 
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ANGLES changes the sine and cosine data from pairs of receivers 
into their phase difference. The sine and cosine of this phase difference 
are first calculated, using the four arguments and the data generated by 
SINCOS. Error 11 occurs if both the calculated sine and cosine are 
zero; otherwise, the sine and cosine of the phase difference are used 
to find the phase difference in the region -tt to it. 

ARRIVL uses the phase differences between pairs of antennas and 
the antenna spacing data to find the angles of arrival . Two error 
conditions exist. If the phase differences between both the azimuth and 
elevation antenna pairs are zero, error 31 is flagged. The subroutine 
calculates the cosine of the elevation angle, and then generates the 
elevation angle itself. If the cosine of the elevation angle is greater 
than unity, error 32 occurs. 

TRGINV takes the sine and cosine of an angle as input, and returns 
the angle itself in radians. The returned angle is in the range -it to tt. 

Table 5.1 illustrates the interaction of the data-processing sub- 
routines . Before a change is made in any of the six programs , the effect 
on all other programs must be considered. 

Four Common Blocks are used by the processing system to store the 
hardware system and echo parameters. Block DIST contains data related 
to the distances between the antennas . This data is used by the program 
ARRIVL. The block holds the reciprocals of the antenna spacings, divided 
by 2tt, since this is the form of the spacing data required by ARRIVL. Also, 
the reciprocal of the distance between the widely -spaced antenna pair, 
half of the reciprocal of this distance, and the distance itself are 



Table 5.1 


Data processing subprogram interaction 


Program Name 

Called By 

Required Subroutines 

SDRIV1 

— 

ANGLES 

ARRIVL 

DOPPLR 

SINCOS 

ANGLES 

SDRIV1 

TRGINV 

ARRIVL 

SDRIV1 

TRGINV 

DOPPLR 

SDRIV1 

— 

SINCOS 

SDRIV1 

TRGINV 

TRGINV 

ANGLES 

ARRIVL 

SINCOS 

— 
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stored for use by ARRIVL. The six values stored in DIST are initialized by 
a block-data program named SPCING. Once initialized, the values are 
not changed. Block ERROR is a single word used for communi eating error 
conditions between the measurement subroutines and SDRIV1. It is initially 
set to zero, indicating no error, by the block-data program ERROR. Block 
TRIG holds the sine and cosine of a, 2a, and 3a, as well as a itself. It 
is used for communication between the programs SINCOS and ANGLES. Block 
AOFA holds the angles of arrival. The block is written by the program 
ARRIVL, and read by SDRIV1. Table 5.2 illustrates how these various 
blocks are used for communication. 

5.2 Simulation Program 

There are two sources of noise in the system. First there is sky 
noise. This noise is Gaussian distributed, and characterized by a 
standard deviation. It is additive with the echo signal at the receiving 
antennas. The second noise source is quantization by the A/D converter. 
This noise is uniformly distributed on the interval [0, 1). Since the 
A/D converter operates on the received signal, quantization occurs on the 
sum of the received echo and the sky noise. The simulation takes into 
account these noise sources by modeling the signal flow through the system. 
The echo signal itself is created first. Then, Gaussian noise is added 
to the signal. Finally, the signal is quantized. Echo signal creation 
requires two sets of parameters: the antenna- spacing data, and the echo 

parameters. The spacing data is simply the distance between the antennas 
in wavelengths. The relative antenna placement shown in Figure 3.1 is 
assumed. The distances between antennas shown in Figure 3.1 are suggested; 
the actual distances may vary by a few percent without any adverse effect. 


Table 5.2 


Data processing subsystem common block usage 
Block Name Initializing Program Written into by 


AOFA 


ARRIVL 


DIST 


SPCING 


ERROR 


ERROR 


ANGLES 

DOPPLER 

ARRIVL 

SINCOS 


TRIG 


SINCOS 


Read by 

ARRIVL 

SDRIV1 

ARRIVL 

SDRIV1 

ANGLES 

SDRIV1 
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The spacing data used in the processing system are the actual measured 
values, rather than the proposed distances. The echo parameters include 
the azimuth angle, elevation angle, velocity, and amplitude of the reflecting 
point. With this data, the phase differences between signals at the 
various antennas and the Doppler shift of the received signals can be 
calculated. The processing system deduces Doppler shift and velocity 
from the phase change per pulse repetition period of the signal at 
any antenna. The simulation assumes a PRF of 400 Hz and a radar wave- 
length of 7.33 m, in order to calculate the phase change due to Doppler 
shift. 

One other factor is required before the echo signal is created: 
the absolute phase at any one of the receivers. Recall that the angles 
of arrival are derived from phase differences between antennas, and Doppler 
shift from phase change with time. But from Section 4, it is apparent 
that algorithm accuracies are functions of the absolute phase. The 
simulation takes as an input parameter the number of trials to be performed, 
as well as the echo characteristics. Starting with an absolute phase of 
zero at antenna one for the first trial, the absolute phase for each 
successive trial is that of the previous trial incremented by 2-n/Mj where 
M is the total number of trials. The processing system results for each 
trial are compiled; after the final trial, the average and standard 
deviation for the azimuth angle, elevation angle, and velocity returned 
by the processing system are output. In this way, the effect of absolute 
phase on measurement accuracy is included in the overall accuracy 


determination. 
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Six quantities specify the characteristics of an echo. These 
quantities are input from the keyboard, and used to construct a series 
of reduced data frames. The first quantity input is the echo amplitude. 

This quantity is specified in terms of the A/D converter output scaling, 
and can vary from zero to 511. It controls the effect of quantization 
noise on processing accuracy, in that the created echoes are scaled to 
this amplitude, and then truncated to an integer value. The second quantity 
input is the standard deviation of the Gaussian noise which is to be added 
to the echo signal. This quantity is also specified in terms of A/D 
converter scaling. Signal -to-noise ratio is set through both the amplitude 
and noise specification. For 20 dB S/N, amplitude may be set at 100 and 
noise at 10 units. The third and fourth quantities input are the azimuth 
and elevation angles of arrival in radians; the fifth is the radial 
velocity of the reflecting point towards or away from the receiving array. 
The last quantity input is the number of trials to be performed with 
the same echo parameters. After the last quantity is input, the simulation 
and processing begins. The specified number of trials is made, each 
with different initial phase, and different random values of Gaussian 
additive noise. The results of the test are printed, and the simulation 
program is ready for a new test specification. A listing of the simulation 
program is given in Appendix V. A sample run is shown in Table 5.3. 

5.3 Simulation Results 

The simulation testing is actually two independent test series. 

The first series verifies the correct operation of the processing system 
under no-noise conditions , and determines the range of values of the echo 
parameters for which the processing system produces the required results. 
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Table 5.3 

Sample runs of the TRY4 simulation program 


************************************* 
AMPLITUDE? 1 TO 5 00) 

>im 

NOISE STAND .DEV. 

>2 

ELEVATION ANGLE (PAD) 

>.7854 

ANGLE WITH DUE NORTH (RAD) 

> .2 

RADIAL VELOCITY (MS) 

>20 

NO. 0F ITERATIONS 
>2 9 

AVERAGE STANDARD DEV 

AZIMUTH 0.2002 0.0085 

ELEVATION 0.7g53 0.0021 

VELOCITY 20,5373 2.3077 

************************************* 
AMPLITUDE (I TO 500) 

>100 

NOISE STAND .DEV. 

>2 

ELEVATION ANGLE (RAD) 

> . 7854 

ANGLE WITH DUE NORTH (RAD) 

>.2 

RADIAL VELOCITY (M/S) 

>50\N2 0 

NO. 0F ITERATIONS 
>29 

AVERAGE STANDARD DEV 

AZIMUTH 0.2 003 0.0073 

ELEVATION 0.7853 0,0018 

VELOCITY 20.2524 2.2158 
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The second series determines the measurement accuracies of the processing 
system when additive Gaussian noise is present, and finds the minimum 

signal -to-noise ratio for which the required measurement accuracies are 
achieved. 

5.3.1 Validity test. In order to test the validity of the processing 
system, a parameter list is given to the simulation program, which then 
produces a reduced frame. This reduced frame is passed to the processing 
system, which attempts to recover the input parameters. For the validity 
test series, zero Gaussian additive noise is used; however, quantization 
noise is present. Therefore, the parameters input are never recovered 
exactly. The criteria used for identifying a processing system failure 
is the following: if any of the recovered parameters vary by more than 

just a few percent from the input parameters, a failure has occurred. 

Using this definition, the input parameters are varied throughout their 
expected range, and the validity region is delineated by the input para- 
meter values for which processing system errors occur. In mapping this 
region, a standard input parameter list is vised. Each parameter in the 
list is varied independently until a processing system error occurs. 

The points of error are noted, and then that parameter is returned to 
its standard value, and another parameter is varied. This process 
is continued until the validity region is completely determined. The 
standard parameter list is as follows: 

Azimuth Angle : 11° 

Elevation Angle: 45° 

Velocity : 50 msec 

Amplitude 


100 
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These parameters are selected in view of the transmitting and receiving 
antenna patterns from Section 3, and the meteor-trail characteristics 
from Section 2. They represent typical operating parameters, rather 
than the optimal conditions for the processing system. The amplitude of 100 
sets the quantization noise at less than one part in 100, as described in 
preceding paragraphs . 

Test results show the validity region for elevation angle to be 0 to 
90°. The error in measurement of this quantity is relatively constant from 
0 to 70°, but rises sharply above 70°. The required validity region is 
0 to 90°, so that satisfactory performance is achieved. For azimuth angle, 
the validity region is -45° to 45° . According to the specification in 
Section 2, the required region is -90° to 90°. The equations in Section 
4.4 indicate that the validity region cannot be increased without decreasing 
the azimuth antenna pair spacing from its current dimension of one radar 
wavelength. The error in azimuth angle measurement is nearly constant 
across the validity region. The velocity measurement is valid from -300 
to 300 msec . The measurement error is constant for wind magnitudes 
under 100 msec , but grows for wind magnitudes greater than 100 msec 
This is the expected result, as indicated in Section 4.2. The processing 
system satisfies all of the measurement region requirements, except for 
azimuth. A discussion of methods for increasing the azimuth validity 
region is given in Section 6. 

5,3.2 Noise test . In order to find the minimum signal-to-noise ratio 
(S/N) for which the required measurement accuracies are achieved, a 
series of tests is conducted. Again, a standard parameter list is used, with 
the following values chosen; 



Azimuth angle 
Elevation angle 
Amplitude 


80 


11 ° 

45° 

100 ° 


These parameters are held constant throughout the test series. The 
velocity and Gaussian -additive noise are varied, and the test results 
tabulated. In a preliminary test series, measurement accuracy was 
found to depend on velocity magnitude but not direction. The values 
of velocity used in the test series are, therefore, chosen as 0, 10, 20, 
50, 100, and 200 msec ^ . The Gaus sian-additive noise values used in the 
test series are 1, 2, 5, 10, and 20 units. Since the amplitude is set 
at 100 units, the corresponding S/N are 40, 34, 26, 20, and 14 dB. 

With 6 values for velocity and 5 for additive noise, a total of 30 tests 
are made. Each test consists of 29 trials. The average and standard 
deviation of azimuth, elevation, and velocity for each set of trials 
are the outputs of each test. 

For all tests, the averages agree closely with the input parameters. 
The standard deviations of the output quantities increase with increasing 
additive noise. Figures 5.1, 5.2, and 5.3 show the standard deviations 
of the velocity, elevation, and azimuth measurements as a function of 
additive noise for various velocities. Note that the curves on all three 
figures are indistinct for velocities less than or equal to 100 msec 1 . 
This supports the previous conclusion that velocities below 100 msec 
do not affect measurement accuracy. The curves on all three figures 
show a marked increase in error for a velocity of 200 msec \ however. 
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Figure 5.3, Elevation standard deviation versus additive noise. 
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The minimum S/N which gives the required measurement accuracies is 
found by considering the 100 msec * curve on each figure. 

From Section 2, the velocity accuracy requirement is 5 msec Figure 
5.1 indicates that a S/N of 30 dB gives this result. The required elevation 
accuracy is .003 rad; from Figure 5.2, 33 dB gives this result. Finally, 
the azimuth accuracy requirement is ,1 rad, which is achieved for a S/N 
of 17 dB. Obviously, the elevation angle is the most critical. 

From Section 1, the highest the S/N can ever be is 40 dB. The 33-dB 
requirement is under this 40-dB limit, and is, therefore, achievable. With 
33-dB S/N, velocity accuracy is 4 msec -1 and azimuth accuracy .01 rad. 

Then all requirements are met by a signal-to-noise ratio of 33 dB. 
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6. RESULTS AND SUGGESTIONS FOR FUTURE WORK 
In preceding chapters, the hardware and software subsystems which 
make up the meteor-radar system are described. Two points are addressed 
in the following subsections. First, what combination of pulse width, 

PRF, and receiver bandwidth is optimal for the total system. Second, 
how can the azimuth validity region be increased. In the last subsection, 
conclusions and suggestions for future work are presented. 

6.1 Pulse Width , Bandwidth, and PRF Tradeoffs 

The system hardware and measurement algorithms place constraints on 
three system parameters: transmitted pulse width, PRF, and receiver 

bandwidth. These constraints are considered individually in previous 
sections. Here, all constraints are brought together, tradeoffs are 
evaluated, and the final parameter selection is made. 

The transmitter is rated for a nominal average power of 20 kw, and 
a maximum average power of 40 kw. The instantaneous power is fixed at 
5 Mw. The effect of fixing the average and peak powers is to impose an 
inverse relation between pulse width and PRF. The products of PRF and 
3-dB pulse width are constrained to .004 for nominal power operation, 
and .008 for maximum power operation. 

Two range algorithms are given in Section 4.1. One uses a ratio method, 
the other a parabola method. The parabola method requires a characteristic 
pulse at the receiver output wider than 30 usee at 20-dB points (as 
opposed to the 20-usec requirement of the ratio method) but gives a more 
accurate range determination than the ratio method, and is used in the 
system. For a Gaussian pulse, 30-usec wide at 20-dB points, the 3-dB width 
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is 12 ysec. With multiple targets, however, the narrower the characteristic 
pulse, the better the resolving power of the radar system. The character- 
istic pulse width is then 12 ysec, and no wider, so that the range 
algorithm requirement is satisfied, and the resolving power is not 
compromised. 

In passing through a finite-bandwidth receiver, the shape of an input 
pulse is modified according to the relations given in Section 4.1. 

Relative to an infinite-bandwidth receiver, the output pulse length is 
increased, and the amplitude decreased. The range- algorithm pulse-width 
constraint applies to the receiver output, rather than the input. In 
order to minimize the signal amplitude loss in the receiver, a bandwidth 
which decreases the signal by no more than 5 percent is needed. Solving 
equations (4.3) and (4.4) (using a characteristic pulse width of 12 ysec) 
the required receiver bandwidth and transmitted pulse width are 88 kHz and 
11.4 ysec. 

The Doppler algorithm gives good results for a PRF in the range 400-1600 
Hz, with optimal results at 700 Hz. From the system interconnection plan 
in Section 3, however, the highest usable PRF is 625 Hz; any higher PRF causes 
range ambiguities. Then the acceptable range for the PRF is 400-625 Hz. 

Using 400 Hz and 11.4 ysec as the PRF and transmitted pulse width, 
the average transmitted power is 23 kw, very nearly the nominal power 
rating. Increasing the PRF to 625 Hz brings the average power to 36 kw. 

A PRF, pulse width and bandwidth of 400 Hz, 11.4 ysec, and 88 KHz 
are reasonable parameters to use initially. Of course, the assumptions 
of Gaussian receiver bandpass and transmitted pulse shape are only 
approximations. The transmitted pulse width may require some stretching 
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or shrinking in order to get a receiver output pulse 30 ysec wide at 20- dB 
points. Once the transmitter is completely operational, the PRF can be 
increased to 625 Hz. At that point, increasing the PRF or pulse width 
will degrade system performance by introducing range ambiguities and 
increasing the likelihood of multiple-echo overlap. 

6,2 Azimuth Validity Region 

With the antenna arrangement of Figure 3.1, the azimuth validity region 
is from -45° to 45° relative to due north. The specifications from Section 
2 call for the region to extend from -90° to 90°. Although few echoes 
are expected from outside the main transmitter beam, which extends from 
-26° to 26°, some means for identifying out-of-beam echoes is needed. 

Several methods to achieve this identification are given below. 

If the azimuth antenna spacing is decreased to X/2, the validity region 
becomes -90° to 90°. This straightforward solution is difficult to 
implement, however, because of mechanical and electrical interference 
problems with such a closely spaced antenna pair. Even the current IX 
spacing is perhaps closer than desirable. This solution therefore cannot 
be used. 

Another available solution uses the high directivity of the Yagi 
antenna azimuth pattern to identify out-of-beam echoes. The azimuth 
antenna pair phase-center spacing is maintained, but the Yagi antennas 
which form the pair are skewed from due north by a half-beamwidth, one 
to the left and the other to the right. Then echoes which are in the 
main beam of both antennas are from the desired region. Echoes outside 
this region are identified by the signal amplitude differential between 
the left-pointing and right-pointing antennas. The main difficulty with 


this approach is in locating the phase center of the Yagi. Very little 
published information is available on this topic. Another problem is 
the gain-match required of the two receivers connected to the azimuth pair 
antennas. Although difficulties arise in implementing this solution, it 
is a practical technique for recognizing out-of-beam echoes. 

A third solution adds a dipole antenna to the receiving array. Echoes 
out-side the main beam are recognized by the relative strengths of the 
signals from a Yagi and dipole. Much stronger signals are available at 
the Yagi than the dipole for echoes in the main Yagi beam. Outside 
the beam the signals from a Yagi and dipole are nearly the same. This 
solution is essentially similar to the previous one, but avoids the 
Yagi-phase- center problem. However, extra multiplexer channels are needed, 
or the multiplexing sequence must be modified, in order to add the additional 
information from the dipole to the frame. 

The azimuth ambiguity problem indicates that more work is needed in 
the area of the receiving array design; perhaps a completely new antenna 
arrangement is needed. The best approach to the problem of antenna array 
design at this point is to implement a portion of the receiving array, 
and determine the extent of array modification needed to satisfy the 
system requirements. 

6.3 Suggestions for Future Work and Conclusions 

In previous Sections, the hardware and software subsystems required 
to build a meteor-wind system for observing short-period dynamics are 
described. The system is built around a medium-scale computer; all meteor- 
echo recognition and parameter measurement algorithms are incorporated in 
software. The computer eliminates many complicated hardware components. 
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which results in an increased mean-time between failures. All meteor-echo 
parameters are measured to the required accuracy for observing gravity 
waves; in particular, a height accuracy of ±1 km is achieved. The system 
uses a high-power transmitter and a high-speed collection system to achieve 
temporal as well as spatial resolution. 

The problem of azimuth-angle ambiguities for echos outside the main 
lobe of the transmitting antenna requires more study. In fact, the entire 
receiving-antenna array should be reevaluated, with the intent of finding 
an arrangement of antennas which has no azimuth or elevation angle 
ambiguities. With such an arrangement, an all-sky system is available 
by simply replacing the high-directivity transmitting and receiving 
antennas with dipoles. The software system is modular, so that a different 
receiving array requires only minor modifications to the system programs. 

System calibration is also an important, unresolved problem. Such 
things as cable delays and receiver signal delays must be calibrated out of 
the system. One method for calibration uses the known path of satellites. 

The radar system tracks the satellite, and the known path of the satellite 
is compared to the path given by the radar system. Another approach uses 
a low-power transmitter suspended from a balloon. The balloon is released 
north of the receiving array, so that it rises through the antenna boresight. 
The transmitter is tracked by the receiving system, as well as by a 
theodolite. The optical measurement is used to calibrate the radar 
measurement. Both these methods are reasonable approaches to the initial 
calibration problem; however, they are not suitable for use on a routine 


basis . 


Although 9 channels of useful data are available at the receivers 
during echo reception, only one channel can be digitized at a time. 

A method for sampling all channels simultaneously during echo reception 
would simplify the data reduction subsystems, and also increase the 
useful information gotten from each meteor trail. With the addition of 
some hardware, the radar system can operate in this mode. The video 
channel is digitized continuously, until a meteor-trail echo is detected. 

Then, under program control, the digital input/output interface to the 
computer is used to output the time in psec after pulse transmission when 
the next echo is expected. A counter, driven by a 1 MHz clock, is reset 
to zero as the next pulse is transmitted. A comparator senses when the 
count equals the expected delay in psec; sample-and-hold circuits connected 
to each of the eight phase -comparator channels are strobed at that instant. 

The A/D converter is then used to serially input the 8 samples. The total 
hardware required for this mode consists of a counter, comparator, and 8 
sample-and-hold circuits. The 1 MHz clock signal is available from the 
A/D converter. Since the required hardware is mainly digital integrated 
circuits, the physical volume and interconnection complexity are low. The 
addition of this small amount of hardware replaces a large amount of 
software, and results in more data collected per transmitted pulse. Modifying 
the system in this way may prove useful, once the basic system is 
operational. 

Ground-based measurements of short-period wind variabilities in 
the meteor region are essential to the investigation of atmospheric 
dynamics. The system described in previous chapters should give a higher 
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echo rate and better resolution than any meteor-wind system currently 
available, and contribute significantly to the understanding of atmospheric 


winds . 
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APPENDIX I 

PRINCIPLE OF DOPPLER DETERMINATION 
Measure the quantities M^, defined as 


Afj = 4 cos (6 + toi ) 


M 2 - A sin(0 + 

= A cos(0 + utj) 

0 = phase angle. 

til = Doppler frequency. 

A = Scale factor corresponding to the amplitude of the echo. 

The three measurements are taken sequentially, with the time between 
measurements given by the pulse repetition period At; therefore 


t 2 " *1 " *3 ' t 2 = At 


The three measurements are then rewritten 


Let 


Then 


Also 


~ A cos (6 + tat 2 - wA t) 
M 2 = A sin (6 + <at-^ 

= A cos (6 + tat 2 + mAt) 
0 = 0 + 10^2 

a = taAt 


«* A cos (6 - ot) 

M 2 - A sing 
M j = A cos (0 + a) 

- A cosg cosa + A sing sina 


= A cosg cosa - A sing sina 



Combine the three measurements as follows: 


2 A sing sina 
2 2 A sinB 

* sina 

= sin coAt 


a) A t 


sin 


-1 


W 1 - 


2 M„ 


_1 

At 


sin 


-1 



(A) 
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APPENDIX II 

QUANTIZATION ERROR IN DOPPLER ANGLE DETERMINATION 

Let M 1 + ^ = A cos(B - a ) 

A ?2 + = A sinfi 

- A cos CP + a) 

i = 1, 3 are truncation errors due to A/D conversion. 


0 < e . < 1 
M - M 

z ' -- -nq for ^2 * 0 


2 A sin$ sina - + e 



Since / 0, - — + 2 -— ^ 1. Then the expression for Z' can be expanded 


+ higher order terms 
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The error in Z' is defined as 


A Z' = Z' - sina 


Then if values for e , e 2 , and are chosen such that the error is 
maximized, an upper bound on the error in Z* is found. 


I AZ 1 


1/2 + sin a 


Mn 
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APPENDIX III 

THE PESSIMUM CONDITION FOR DOPPLER DETERMINATION 


Let 


Then 


^2 = A sin(0 + 

- A cos ( 6 + oitj) 

<j> = 0 + 

a - u>(t 3 - t^) 

= A sin $ 

Afj - A cosO + a) 


Given a, find the lowest magnitude that the higher of the two 
measurements M 2 and can reach for arbitrary values of <fr. The pessimum 
condition is that \M 2 \ = |^ 3 I- 


Assume 


\M 2 \ = |Af 3 | 

| A sin<{>| = \A cos(<t> + a) | 

|sini})| = |sin(± 21 + 4> + a] | n =1, 3, 5, . . . 
±4> = (± 21 + + a) 


There are four solutions : 


( 1 ) 

(2D 


, . MT 

(|) = ±^— +4> + a 
— nir 

“ = + “ 


-4> 

♦ 


^ nir 

± 2 ~ * <t> + a 


— K7T 

+ “ 


a 

2 
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(3) 

(4) 



-♦ 


+ 



a 


+ TL 

~ 2 


j^ZTT 

Cases Cl) and (4) , where a = ± j- , imply that the minimum of M 2 or 
is zero. Cases (2) and (3) are the required solutions for -it/ 2 < a < tt/ 2. 
Cases (2) and (3) can be combined 



n = 1, 3, 


Then the pessimum condition gives that 


(a) 






M 2 

= 

A sin 

7T | a | 

4 ' ' 2 ' 

K J 


or 


00 







_ 

A sin 

|l- ItI 

2 


4 1 2 1 




w , 


Consider case (a) only. Case (a) is the pessimum condition for 
-u/2 < a < ir/2: case (b) is the condition for it < |ct| < n/2, which 
is not of interest. Then the lowest magnitude that the higher of \Mj\ 
or 1*3 1 can have is 


Max 



A sin 


I l2l 

4 2 
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APPENDIX IV 

INTERFEROMETER PHASE DIFFERENCE DETERMINATION 


Measure at time 

ty = Aj sinC^j + wtj) 

t 2 : m 2 = A i cos ^i + 2 5 


from Antenna #1 


= A 2 sin((J> 2 + w* 3 ) 
*4 : = A 2 cos ^2 + W V 


from Antenna #2 


Let 


tut, = 0, 


tut- = 6 


B 


tut- = 0, 


tut. = 0, 


Then = A^ sinf^ + 0^) 

* A, sint}). cos0. + A. cost)), sin0. 

11 All A 

M 2 = A x cos (4^ + 0 B ) 

= costfrj cose^ - A^ sintfij sine^ 

= A 2 sin(<|> 2 + 0^) 

* A 2 simj> 2 cose^ + A 2 cos(J> 2 sine^, 

A/ 4 = A 2 cos(<t > 2 + Q D ) 

= A 2 cos<f> 2 cosQ^ - A 2 sin4> 2 sin8^ 



M cos0 D - M sine. 

1 d 2. A 

np— 

cose + M sin0„ 

2 A JL D 

cos (6^ - e B ) 

M- cosG n - M. sine„ 

3 D 4 C 

coste^ - e^) 

M 4 cosG^, + Af 3 sinS^ 
coste^, - y 

Let X = ^4^2 cos^ - 4> 2 ) cosfe^ - 6 fi ) cos (0^ - 0^) 

Y = i4jA 2 sin((|)j - ij> 2 ) cos (8^ - 0 fi ] cosCe^ - 6^) 

I = [A^ cos^j A 2 cos<j> 2 + A 1 sin<j> 1 A 2 sin^] cos(0^ - e g ) cosfG^ - 0^) 
= sinG^ + M 2 cos0^][Af 3 sinG^ + cosG^,] 

+ [ M ^ cosG^ - M 2 sin0^][A/ 3 cosG^ - M 4 sinG^] 

X = M ] M 5 cos(0 g - e D ) + M 1 M 4 sin (0^ - e^) - sin(G^ - d D ) 

+ m 2 m 4 cos(e^ - e^p 

Y = [A^ sin<t> 1 A 7 cos 0 ^ cos^ A 2 sincj> 2 ] cosC©^ - 0^) cosCG^, - 0 ) 

= [M 1 cos0 B - M 2 sinG^] [M^ sinG^ + M 4 cose^,] 

- [M 5 cos& d - M 4 sinG^] [M^ sine^ + M 2 cosG^] 

Y = M y M 4 cos (0£ - B c ) - M 2 M z cos(G^ - G D ) - M ] M l sin(G B - 9^) 

- M 2 M 4 sin(0 4 - G^) 

If 6 D - B C ’ e c - °B = e B - 9 A ' “ 4t 
and wAt = a 


A x sin<t> 1 = 
A ^ cos^ = 
A 2 sin<^ 2 = 
A 2 cos(j) 2 = 
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X - (A/jA/j + A^A/^) cos ^ a " sino + sin3a 

Y = (MjA/j + A^A/^) sin2a + A/^A^ cosa - cos3ct 


From the definitions of X and Y, 


AT 2 + Y 2 ■ A* 4 2 cos 2 ^ - 0 fl ) cos 2 (0 c - 0 D ) 


sinC+j - <fr 2 ) = 


/ A + Y* 


cos (4^ - 4 > 2 ) 


/Y 2 + Y 2 


Note that these ratios are singular only when cosot = 0. 



o o O O 
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APPENDIX V 

METEOR RADAR PROGRAMS 


C. BACKOF 
C EDIT #14 MO D2 

C THY 4 IS A PROiRAM TO SIMULATE THE DATA FROM THE METEOR 

C TRAIL SEAHCH AND STORE PR03RAM. IT TAKES AS INPUT DATA 
C THE SYSTEM DESCRIPTION AND THE COORDINATES AND DELOCI TY 
C OF REFLECTION POINTS. ALSO THE NUMBER OF ITERATIONS IS 

G INPUT* THE INITIAL PHASE IS ROTATED 3Y 2PI/N RADIANS FOR 

C EACH TEST* WHERE N IS THE NUMBER OF ITERATIONS. THE RES- 

C ULTS OF EACH TEST ARE ADERA3ED* AND ONLY THE A DERATE AND 

C STANDARD DEVIATION ARE PRINTED. THIS TEST DOES NOT IN- 
C CLUDE THE RANT E ALJORITHM. 

C 

APRIL 19 1 A 73 
GUT SYSTEM DESCRIPTION 

DIMENSION NCI 6), A N( I«),K(32) 

yPlTECSl) 

FOPMATC !H1 # 22X,2*HT P Y OUT PRA.MF GENERATOR /// 

11 H ,P- 7 X, 1 PHSYSTEM DESCPIPTIOM //1?H AZIMUTH SPACI*<0: ) 
CALL FFIP ( 1 ,AZSP > 

2 FORMAT C FI 0.3) 

’•'P I TF (33) 

3 FOP MAT ( //S4H FINE ELEVATION SPACING; ) 

CALL PFIRCl ,FELSP) 

VP I TEC 3, 4) 

4 FORMAT C//P3H COARSE ELEVATION SPACING. ) 

CALL FFIP Cl, CELS?) 

I NT- P2 P7 

C 

c ASK TO 55 PAPTICULA°S OF EACH METEOR 

r 

5 U P ITE (3,3) 

< FORMAT (3 PK ************************************* 

1/21H* AMPLITUDEC ! TO OOR). ) 

CALL FFIP (I, A) 

VP I TEC 3, 7) 

7 FORMAT (17H ”01 SF STAND .DEV. ) 

CALL FF IP C l , SD) 

' ,? P ITE C 3, P) 

p FORMAT (22H ELEVATION ANGLE (PAD) ) 

CALL FFIP(1,FL) 

VRITEC 3,9) 

0 FORMAT (2«H ANGLE WITH DUE NOPTH (R9D) : ) 

CALL FFIR ( 1 ,A Z) 

VP I TEC 3, l?) 

10 FORMAT (23 H RADIAL VELOCITY CM/S) ) 

SUMAZ=0. 

SUMEL=7I. 

SUMVLrO. 



SSOAZr?, 

SSOELrO. 

SSQVL=«. 

CALL FFIRCi.V) 

'•’PI TEC 6, 1?) 

13 FORMAT C IBM NO. OF ITERATIONS ) 

C 

CALCULATIONS 

C 

WDT = V/233.2036 
CALL FFIRC1,«U) 

IQ=nu 

CSEL= COS (EL) 

FELD =- 6.2 R32* EEL SP*CSEL* COSCAZ) 
CELD=6.2B32*CELSP*CSEL* COSCAZ) 
AZD=-S.2P32*AZSP*CSEL* SINCAZ) 

HO 133 11=1,10 
PH =6 .2 B32/QU*EL0AT (II) 
no 15 1=2,5 
15 A NC I ) =PH 

C A HOLES AT ANTENNAS 

DO 17 1=6,9 
1? AN(I)=PH+FELD 

A NCI) =3. 

A N(!3)=0. 

ANC11)=PH+CELD 
A NC 12) =A NC 1 1 ) 

A N< 13) :PH 
A MO A) = D H 
A M(15)=PH+AZD 
AN( 1 5) =A N(1 5) 

r* 

C PUT IN DOPPLER FREQUENCY CHANGE 
DO 20 1=1, 1$ 

A NCI) -A MCD + WT* FLOAT Cl) 

20 ANCI)=AM0D(ANCI),6.2B32) 

00 30 I =2 9 2 

30 NCI ) = A* SI NCA NC I) HGAUSSCI NT)*SD 

DO 40 I =3, 9,2 

40 NCI) =A *COS C A NC I ) )+ GAUSS C I NT) *SD 

DO 50 1=11, 1C, 2 

50 NCI)=A* SI NCANCI)>+ GAUSS Cl NT)* SD 

DO S3 1=12, 16,2 

60 NCI)=A*COS(AN(I)) + GAUSSCINT)*SD 

MC l)sSS>D/3.+ .5 
NC 13)2 NCI) 

CALL SDRIVl ( N,AZ ,EL, VEL) 

SUMAZ = SUMAZ+AZ 
SUMELr SUMEL+EL 
SUMVL = SUMVL+ VEL 
SSQAZ=SSQAZ+AZ*AZ 



SEOEL = SSOEL+ SL* SL 
SSOVLr SBQVL+VFL* '/EL 
AVAZrSUMAZ/HU 
AVEL=SlJMEL/OU 
AVVL = SUMVL/Q'J 

SD AZ = SOFT ( (SSfiAZ -SUMAZ* SUMAZ /nil) /(OU - 1 •) ) 
SDELr SOFT C (SSOEL-SU MEL* SUMEL/OU) /(Of-l .)) 
SDVLrSCTPTC (SSOVL-SUMVL* SUMVL /nu> /COU- 1 . > > 
VF I T E < fi . , 2 0 3 ) A' VQZ , S 9 A Z , A ''EL , SO EL » A VVL , S D VL 
FORMAT C12X # 7HAVFRA0E f Sy t 1P.HSTA NOARO OEV/ 

1 1PIH AZIMUTH ,F9.4,7X,FP.4/ 

2 lflH ELEV?3TI0M # F9.4 t 7X t F«?.4/ 

• 3 I AH VELOCITY , F9 .4 1 7X, FP .4) 

OOTO *5 
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BLOCK DATA 
C. BACKOF 
C ED IT #1 

C ERROR INITIALIZES IERR TO ZERO. . e THAT’S ALL! 
COMMON /ERROR/ IERR 
DATA IERF/fl/ 

F ND 


SLOCK DATA 
C. BACKOF 
C EDIT #1 
C .JULY 19*1973 

C SRC I NO IS THE PROGRAM WITH THE ANTENNA SPACING DATA 
C 

COMMON /DIST/ AZSP,CELSP,FELSP,R2PI ,RPI ,RR2PI 

r, 

C THESE ARE THE PECIPROCLS OF ?. PI* ANTENNA 5PACINGS 

r* 

DATA AZSP.CELSP.FELSP/. 1591549,. 1591549, 7, 957745E-3/ 

C 

C THESE ARE 1 /FELSP , I /2 /FELSP , FELSP 
C 

DATA R2PI f RPI ,RP2PI /. 05 f . 025-, P.0 . / 

END 


SUBROUTINE ANGLES ( Nl , N2 , N 3, N4, ANGLE) 

C. BACKOF 9 MARCH 73 
C EDIT #11 

C SUBROUTINE ANGLES TAKES INFO FROM THE SINE AND COSINE 
CHANNELS OF TWO RECEIVERS AND GIVES THE SINE AND COSINE OF THE 
C PHASE DIFFERENCE 

COMMON /TRIG/ SI N1 ,S I N2 , SI N3,C0SI ,C 0S2 ,C0S3, ALPHA 
COMMON /ERROR/ IERR 
C 

CALCULATE FLOATED VALUES 
C 

A 1 = FLOAT ( N1 > 

A2rFL0ATC^) 

A3 = FLOAT ( N3) 

A4r FLOAT ( N4) 

C ' 

CALCULATE INTERMEDIATE OUANTITIFS 

C 

01 =A 1*A4 
02= A !*A3+A2*A4 
03 =A2*A3 
C 

CALCULATE UNNORMALIZED SI NE A ND COSINE 

C 

C0S=O2*C0S2-Qi*SI N1+Q3*SI N3 
SI N=02*SI N2+QJ *C0S1 -Q3*C0S3 



c 

C TDST FOR ERROR CONDITION 

Q 

I F( SIM .FO. 0. .AMD. COS . FO . 0. ) GOTO 1 
C 

CALCULATD ANGLE 
C 

A NGLE=TRG I NV( SIM, COS) 

I F( ANGLE ,LT . 0. ) A.NGLE=A NGLE+F.RP32 
C 

c FXIT TO CALLING ROUTINE 

C V 

RETURN 

C 

C RPROR FXIT 
C 

1 I ERR = 11 

RETURN 
END 


SUBROUTINE ARRIVL 
C. BACKOF ?, a MARCH 73 
C EDIT *11 

C SUB. APRIVL TAKES PHASE DIFFS AMD ANTENNA SPACINOS, AND 
C RETURNS ANGLES OF ARRIVAL 
C 

C SET UP ANTENNA SE PER AT I ON DATA 
C 

COMMON /DIST/DAZ t DCEL t DFEL,P2PI ,P D I,PPP°I 
COMMON /A OF A /OAF ,OCEL , OPEL, AZA NO, EL A NG " 

COMMON /ERROR/ IERR 
C 

C ^QPMALIZE PHASE 
PHAZ=OAZ*DAZ 

IFC PHAZ .GT. .5 ) PHAZ=PHAZ-!. 

PHCEL=OCEL*DCEL 
PHFEL=OFEL*DFEL 

COARSE AND FINE ELEVATION 

PHEL=PHFEL+R2?I*AI NT (PHCEL* RR?. p I ) 


ERR =PHEL-PHCEL 
I F ( ERR ,LT. -RPI ) 

I F( ERR .GT. RPI ) 

C 

CALC SUM OF PHASES SQUARED 
C 

CELSQ=PHEL*PHEL+?HAZ*PHAZ 


?HEL=PHEL+R2PI 

PHEL=PHEL-RRPI 


COMBI ME 


CORRECT 

C 
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C 

CALC TRIG 

o 

0 

IFC CEL SO ,EQ . .Cl) CO TO 1 
CEL = SORT C CEL SQ) 

IFC CEl.pQ .GT. 1. ) CO TO 2 
SEL= SORT ( I ,-CELSQ) 

C 

CHANGE TO A MOLES 
C 

AZA NG=TRGI NVCPHAZ ,?HEL) 

EL A MG = TR G I N V C SEL , CEL ) 

C 

C EXIT 
C 

RETURN 

C 

C ERR OR EXI T 
C 

1 IEO p r31 
RETURN 

2 I EPF =3?. 

RETURN 

END 


FU NOT 1 0 N DOPRLR C Ml # N 2, NX , N A ) 

C, RACK OF 9 MARCH 73 
C EDIT #F 

C FUNCTION DOPRLR TAKFS 4 SEOUE NTIAL VALUES AND CALCULATES 
0 SIN OF THE DOPPLER ANGLE ALPHA 

r* 

C DEFINE THE TRANSFORMATION 
C 

COMMON /ERROR/ IF.PR 

CRUMCHC ,J,K)=FL0ATCK-I )/FL0AT(J+J) 

C 

COMPARE MAGNITUDE OF ® A N3 

C IF M2 & NX ROTH ZFRO, ERROR EXIT WITH IERR=2!,FLSF 
CHOOSE THE REITER METHOD DEPENDING ON THE GREATER MAGNITUDE 
C 

MAG2 = IA3S ( *£) 

MAGX=IARS ( M3) 

IFC MAG2+MAG3 *EQ. 0 ) GO TO I 
C 

CALCULATE DOPPLR 
C 

IFC MAG2 .LT. MAG3 ) DOPPLR rCRU NCH ( M4, N3 , N2) 

IFC MAG2 .GE. MAG3 ) DOPPLR =CRU NCH C Ml , M2, N3) 

C 

C EXIT TO CALLING PROGRAM 
C 


RETURN 



O O O O O o Ci o o o o o o o o o o o o o o o o o — o o o 
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EPP OR EXIT 

IFP!® =2 1 

RETURN 

FMD 

. BACKOF 
ffr p2 f 74 

THIS FUNCTION CALCULATES TAU FOR THE RA NOE ALGORITHM 
USING THE PARABOLA METHOD 


THREE ARGUAME MTS APE REOUIRED : THE HIGHEST SAMPLE, AND ITS 

TVO HE I G HP OPS 

FUNCTION ITAUC Nl , N2, N3 ) 

I TAUrl FI X (5 .* FLOAT ( Nl -NS ) /FLOAT C Nl + M3 - M2 - M2 ) + . 5 ) 

RETURN 
E HD 


SU3P0UTI ME 9)RI VI ( N f A 1 ,A2,A3) 

. BACKOF 20 APRIL 73 
EDIT ?il MOD 2 

SDR I VI IS THE SUBROUTINE DRIVER. IT FORMATS DATA AND 
CALLS THE PROCESSING SUBROUTINES IN THE PROPER ORDER , 

GET SYSTEM DESCRIPTION 

DIMENSION M(1S) 

COMMON /AOEA/ AZ ,CEL ,FE1. ,AZANG,FLANG 
COMMON /TRIG/ SI ,S2, S3, CI,C2,C3, ALPHA . 

COMMON /ERROR/ I ERR 

TST DOPPLER 

SDO° 1 =D OPPLR ( N( 2. ) , H ( 3 ) , M ( 4 ) , M ( 5 ) ) 

I F( I ERF .ME. 0 ) GO TO 2D 
SD0P2=D0 DO LR (N(S),N(7),M(P) f M(9>) 

I E( IERR .ME. 0 ) GO TO Z* 

SI =(SD0P1+SD.0P2) /2 . 

TST SI NCOS } GET SINES * COSINES 

CALL SI NCOS 

IF( IERR .ME. 0 ) GO TO 20 
CALCULATE VELOCITY 
C 


VEL =2 3 3 .2 03 <P* AL P HA 
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n 

c GET FINE ELEV PHASE DIFFERENCE 

0 

CALL AMGLES(N(4),N(5),N(0,N(7) f FEL) 

IF( I ERR . NE. H ) GO TO 20 
C 

c GET COARSE EL PHSE DIEF 
C 

CALL AMGLESC N( 1 1 ) , N( 12) , N < 13) , N < U ) , CEL) 
IF ( I ERR . NE. 0 ) GO TO 20 

0 

c AZ ANG PHSE DIFF 

CALL AMGLESC N( 13), N(14),N( 15), N(1S),AZ) 

I F< I ERR .NE. 0 ) GO TO 20 
C 

c GET OF ANGLES OF ARRIVAL 
C 

CALL ARRIVL 

I F( I ERR .NE. 0 ) GO TO 20 
A! =AZANG 
A2=ELANG 
A3 = VEL 
C 

C OUTPUT RESULTS 
C 

RETURN 

9G VRITE(S,11) IERR 

I ERR =0 

11 F0RMATC/3X, BH*** ERROR ,13,3 H*** ) 

RETURN 

END 


C. 

C 


SUBROUTINE SI NCOS 
RACKOF 14 MARCH 73 

SUBROUTINE SI NCOS TAKES SI N OF 
SIN AND COS OF ALPHA, 2*ALPHA , 


ALPHA AS INPUT AND RETURNS 
3*ALPHA ,@ND ALPHA ITSELF 


COMMON /TRIG/ SI N1 ,SI N2 , SI N3 ,C 0S1 ,C 0S2 , C 0S3 , ALPHA 

P 

CALC SIN A COS SQUARED 
C 

COMMON /ERROR/ IERR 
SINSQUSI Nl*? INI 
COSSQI= 1.-SINSQ1 
I F( C0SSQ1 ,LT. (0. ) GO TO ! 


C 

CALC COSINE 
C 


COS 1 = SORT (COSSQ 1 ) 



C' o o ouo- o o o o o o o o o o o o o 


c 

CALC ALL OTHER QUANTITIES PROM ADOVE RELATIONS 
C 

SI N?.:2 ,*SI N1 *C0SI 

COSP. = COSSCl -si MSOI 

SI NSsSI N1*C3.*C0SS0! -SI NSQ1) 

C0S3rC0S!*(C0SSQl-3.*SI NSQI) 

c 

CALCULATE ANCLE 

n 

ALPHA =TRSI NVCSI Ml ,C0S! ) 


EXIT TO CALL I NO RROCPAM 
RETURN 


ERROR EXIT 


I ERR =4 1 

RETURN 

END 


FUNCTION TEC I NV( A t B ) 

RACKOF P.l JULY 73 
VERSION EDIT #2 

FUNCTION TPOINV TAKES 2 VALUES AMD RETURNS THE ARCTAN OF 
THEIR PATIO 

HECK FOR ZERO DENOMINATOR 

IFC R .ED. 0. ) CO TO 1 
TR0INV=ATAN2 (A ,P) 


EXIT 

RETURN 


1 IFC A ,GE. 0. ) TRGINV=l.i™79S3 

IFC A .LT. 0. ) TRGINV=-!.579!79« 
RETURN 
END 
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/ 

/ C. BACKOF 
/ FEB P.7,7 4 
/ 

/ CLOCK STARTS THE FIRST WORD IN THE COMMON 
/ SLOCK ‘STATUS* COUNTING SECONDS 
/ 

/ CALLED FROM FORTRAN BY 
/ ‘CALL CLOCK* 


/ 

.GLOBL 

CLOCX, STATUS 


CLOCK 

0 

DZM* 

STATUS 

/ZERO AT START 


.T I MER 

GO, INC ,5 

/FIRST CALL 


JMP* 

CLOCK 


I NC 

0 


/INTERRUPT SERVICE 


DAO 

SAVE 

/SAVE ACC 


.TIMER 

GO, I NC, 5 

/RESTART TIMER 


ISZ* 

STATUS 

/AND COUNT 


LAC 

SAVE 



.RLXIT 

INC 


SAVE 





.END 




FUNCTION GAUSS (II ) 



C GAUSS GENERATES RAND NUMBS WITH A GAUSSIAN DISTRIBUTION 
C C. BACKOF 
C AUG 25 73 
C 

GAUSS =-S. 

DO 1 1=1,12 

1 GAUSS =GAUSS+R A NDU (II) 

RETURN 

END 


C FUNCTION TO CALCULATE RANDOM NUMBERS 

FUNCTION RANDU(IR) 

IR = IR*G53 

IECIR.LT.0) IR =IR -1 
IR=IR + 257 

IFCIR.LT. PI) IR=IR + 131071 + 1 
RANDUr FLOAT (IR)/131 071.0 
RETURN 
END 



Odd o o o o o o o o 1 — O o o o o O o o o o o o o o o o o o o o 
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RACKOF 
* 7 ? 27,74 

meteor echo search 

FIRST, THE NAME 


SUBROUTINE SPCH 


and preliminaries 

COMMON /STATUS/ ITIM ,ILI N ,1 FRM,I TOP, IBOT 
COMMON /FRAME/ I (3200) 

COMMON /STACK/ -J(220) 

DATA ITH/H-V 

FIGURE OUT WHICH PR a ME TO SEARCH 


I c 

I F ( I FP M ,E0 . 0 ) IS=1*R2 

IE=IS+9R 

TFST FOR BAD FRAME 
(IF FIRST WORD IS -I) 

IF( I (IS-1) .EO. -1) GOTO 4 

SEARCH RIPST RANGE LINE ROP MAX 

MAX -0 
K = I (IS) 

I F( K.LT.MAX .OR. K.GT.M1 ) GOTO 2 

MAX-K 

I?=IS 

IS=I S + l 

I F ( IS .NE. IE ) GOTO t 

MAX IS THE VALUE OF THE MAX, IP ITS POSITION 
OFT CORRESPONDING ELEMENT FROM OTHER RANGE 
LINE; GET RID OF CHANNEL ID 

MAX2 = I (IR+000) + l 14SP0 

AND TEST FOP THRESHOLD 

•I F( MAX.GT.ITH .AND. MAX2 .GT.ITH ) GOTO 3 


THRESHOLD TEST FAILED 


RETUR M 



4> O O O O O O W > O O O O o O Cm o o o 


build reduced frame 

K =1 T0P*22+1 

K IS THE FIRST WORD OF THE TOP LINE ON THE STACK 
IF( JCK> •NE.fl 5 GOTO 4 
STACK FULL 


J(K)=I°+1D0-IE 
JCK+I ) = I (IP-1) 

J(K+2) = M4X 
J(K+3)=I (IP+1) 

J (K+4)=I (IP+K99) 

J (K+5 ) =1 (I D +900) 

J (K+F) =1 (IP+9P1) 

J (K+7)=I (IP + 1 0 a!) 

,] (K+B) =1 (IP+20P) 

J(K+Q)=I (IP+30R) 

J(K+10)=I (IP+40O) 

J CK+I 1)=I (IP+500) 

J (K+12) rl (IP+KflR) 

J (K+13 ) =1 (IP+700) 

J (K+14)=I (IP+EOO) 

J(K+15)=I (IP+IOPIO) 

J (K+1F)=I (IP+llOO) 

J (K+ 1 7) =1 (IP+12^0) 

,J(K+lP)=I (IP+13^0) 

J(K+19)=I (IP+1 4 00) 

J(K+20)rI (I P+1500) 

J(K+21)=ITIM 

REDUCED FRAME IS ON STACK 
UPDATE TOP OF STACK 

ITOP=I TOP+1 

I F( IT OP .CE. 10 ) ITOP=C 

AND LEAVE 

I (IE-99) = 0 
RETURN 
END 

/ 

/C. RACK OF 
/ FEB 2*, 74 
/ 

/THIS IS THE ASSEMBLY PORTION OF 
/THE COLLECTION AND R@NGE SEARCH SUBSYSTEMS 
/TO START, USE THE FORTRAN STATEMENT 
/ 'CALL ADC’ 
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/ 


.OLOSL STATUS, FRAME, ADC 
ADC 0 



LAC* 

(IDS 

/.SETUP ROUTINE 


DAC 

HOLD 



JMS* 

HOLD 

/CO TO IT 


703 70 1 


/SKIP INSTRUCTION 


I NT 


/I NT SVC ROUTINE 

/ 




/ 

INDICATE FIRST 

FRAME 

RAD 

/ 





LAW 

-1 



DAC* 

FRAME 


/ 




/FIRST CALL r 0R 

I MPUT 


/ 





LAV 

-144 

/1 00 (DEC) WORDS 


DAC* 

(20 

/TO WORD COUNT 


LAC 

FRAME 

/FRAME +100 


TAD 

(144 

/IS 2ND LINE 


DAC 

START 

/FIRST RUEE 


TAD 

(-1 



DAC* 

(27 

/TO DATA CHAM A DDR 

/ 




/ 

FIX UP ADDRFS 

SES O r 

STATUS WORDS 

/ 





LAC 

STATUS 



I AC 




DAC 

LINE 



I AC 




DAC 

FPM 


/ 




/ 

START COLLECTION AND 

EXIT 

/ 





703 72 A 


/AD WRITE I MIT 


JM°* 

ADC 



/ 

/ INTERRUPT SERVICE ROUTINE 
/ 

I NT 91 

/ 

/ SAVE REGISTERS 
/ 

DAC SAVAC 

LA CO 

DAO SAVMQ 

OLK 

/ 

/ TEST FOR ERROR 
/ 


703 72 1 
JMP 


OVER 


/SKIP ON AD ERROR 



LAC* 

FPM 


SZA 
. DEC 

LAC 

(1SOO 


.OCT 

TAD 

FRAME 


DAC 

HOLD 


LAW 

-1 

/ERROR RLAC TO 1 ST 

DAC* 

HOLD 

/WORD OF FRAME 

/ 

/FIGURE OUT WHICH CHANNEL 

WAS COLLECTED 

/ 

OVER LAC* 

START 

/OFT FIRST SAMPLE 

CLL 


/SHIFT 14 (DEC) RIGHT 

LRS 

1 S 


IAC 


/+1 

DAC 

CHAN 

/CHANNEL ID 

TAD 

(-2 0 

/EQUAL 1 C(DEC) ? 

SZA 

JM° 

DOWN 

/NO 

/ 

/ FRAME JUST 

FILLED 


/ 

LAC* 

FRM 

/COMPLEMENTT 

SZA! IAC 

/FRAME POINTER 

CLA 

DAC* 

FRM 


DZM 

CHAN 

/BACK TO CHAN 0 

.TIMER 

MNG.6 

/START RANGE SEARCH 

/ 

/ FIGURE THE 

START OF THE 

NEXT BUFFER 

/ 

DOWN LA C* 

FRM 

FRAME POINT 

SZA 

LAC 

(20 

/TIMES IS(DEC) 

ADD 

CHAN 


DAC* 

LINE 

/NEXT LINE 

MUL 

144 


/TIMES 1 00(DFC 

LACQ 

TAD 

FRAME 


DAC 

START 

/START OF BUFFER 

/ 

/ NOW SET UP TRANSFER 


/ 

TAD 

(-1 


DAC* 

(27 

/BUFF ADDR 

LAW 

-144 

/WORD COUNT 

DAC* 

(2S 
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/ 

/RESET 

REGISTERS 



/ 

RAR 




LAC 

SAW 



LMQ 




LAC 

SAVAC 


/ 

/ AND 

EXIT 



/ 

703 724 


/ENABLE COLLECT 


DPR 1 


/LEAVE API 0 


J"!P* 

I NT 

/EXIT 

/ 





/ CALL TO FOR TP AN RANGE! SEARCH 


/ 

PNG 0 

JXS* SPCH 
.PL XI T PNG 

/ 

/ STORAGE AREA 
/ 

LINE P» 

PR M 0 

HOLD A 

SAVAC 0 

SA V7!D 0 

START <? 

C HA N 0 

.END 



